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ABSTRACT 


A  conceptual  space  platform  consists  of  a  payload,  a  power 
conditioning  unit  (PCU) ,  and  two  radiators:  the  main  radiator 
and  a  secondary  radiator.  A  computer  program  was  written  to 
determine  the  required  size  of  the  two  radiators  and  the 
temperatures  of  the  PCU  and  payload  for  a  given  platform  power 
level.  An  iterative  approach  is  necessary  because  the  required 
size  of  the  main  radiator  depends  on  the  size  of  the  secondary 
radiator  and  vice  versa.  Also,  the  temperatures  of  the  payload 
and  PCU  depend  on  the  size  of  the  radiators.  The  program  user 
can  subdivide  the  two  radiators  into  any  number  of  nodes  to 
increase  the  accuracy  of  the  radiant  heat  transfer  solution.  The 
use  of  more  nodes  also  allows  better  prediction  of  the  nonlinear 
temperature  drop  that  occurs  across  the  radiators  as  the  working 
fluid  deposits  the  platform's  waste  heat  in  the  radiator.  View 
factor  expressions  are  automatically  calculated  for  different 
choices  of  the  number  of  nodes .  The  user  can  also  select 
different  separation  distances  between  the  various  platform 
structures.  A  model  is  included  to  couple  the  radiant  and 
conduction  heat  transfer  that  occurs  between  the  payload  and  its 
meteoroid  shell  and  between  the  PCU  and  its  shell.  Also,  the 
program  allows  the  use  of  a  refrigerator  to  cool  the  payload.  If 
a  refrigerator  is  used,  the  program  determines  the  amount  of 
additional  thermal  power  needed  to  run  the  refrigerator.  The 
results  of  parametric  calculations  are  included  to  demonstrate 
the  use  of  the  program. 
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1.0  INTRODUCTION 


Proposed  multimegawatt  (MMW)  space-based  platforms 
containing  payloads  of  electronic  equipment  will  require  large 
amounts  of  electrical  power  for  operation.  The  consumption  of 
electrical  power  by  the  payload  will  result  in  the  generation  of 
waste  heat  which  must  be  rejected  to  space,  via  radiant  heat 
transfer,  to  keep  the  operating  temperatures  of  the  equipment 
within  acceptable  limits.  An  active  cooling  system  can  be  used 
to  cool  the  electronic  equipment  whereby  a  working  fluid  is 
pumped  through  a  heat  exchanger  within  the  payload  and  then  to 
an  external  (secondary)  radiator.  The  size  of  the  required 
secondary  radiator  will  depend  on  several  factors  including  the 
amount  of  waste  heat  to  be  rejected,  the  operation  temperature 
of  the  radiator,  and  the  proximity  of  the  radiator  to  other 
platform  components  such  as  the  payload  and  the  power  supply 
(main)  radiator. 

Determining  the  size  of  the  secondary  radiator  is  not 
straightforward  because  the  presence  of  all  of  the  major 
platform  structures  must  be  considered  when  performing  the 
radiant  heat  transfer  calculations.  Figure  1.1  provides  a 
schematic  diagram  of  the  four  major  structures  comprising  the 
platform.  The  four  structures  are  (1)  the  main  radiator,  (2)  the 
secondary  radiator,  (3)  the  power  conditioning  unit  (PCU)  and 
its  shell,  and  (4)  the  payload  and  its  shell. 

The  main  and  secondary  radiators  are  represented  as  disks 
facing  each  other  and  separated  by  some  specified  distance.  A 
disk  radiating  from  both  sides  represents  the  most  effective  use 
of  heat  transfer  area  for  a  single  radiator.  However,  if  more 
than  one  radiator  along  with  other  relatively  hot  components  are 
included  on  the  platform,  this  configuration  will  not  provide 
the  most  effective  use  of  radiating  area.  Besides  the  relative 
radiator  orientation,  the  optimum  radiator  configuration  depends 
on  other  factors  such  as  platform  stability,  maneuverability, 
survivability,  and  structural  considerations.  (The  design  of 
such  a  platform  will  certainly  be  a  complicated  and  challenging 
endeavor.)  The  platform  configuration  used  for  the  analyses 
described  in  this  report  was  chosen  for  modeling  simplicity  and 
to  demonstrate  that  the  position  of  the  platform  components 
relative  to  each  other  strongly  affects  the  radiant  heat 
transfer  calculations.  These  effects  are  demonstrated  by  varying 
the  component  separation  distances  in  the  radiant  heat  transfer 
calculations. 

The  main  radiator  must  reject  the  waste  heat  associated  with 
the  primary  power  source  of  the  platform.  Most  likely,  a  Brayton 
or  Rankine  thermodynamic  cycle  will  be  used  for  the 
thermal-to-electric  conversion.  The  thermal-to-electric 
conversion  efficiency  will  therefore  be  on  the  order  of  20%  to 
30%  and  most  of  the  primary  power  will  be  converted  to  waste 
heat.  For  a  Brayton  cycle,  the  temperature  at  which  power  is 
rejected  to  space  is  not  constant,  thereby  further  complicating 
the  task  of  determining  the  radiator  sizes. 
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The  secondary  radiator,  as  already  mentioned,  is  required  to 
remove  the  waste  heat  associated  with  the  cooling  °f  the 
electronic  components  comprising  the  PCU  and  payload.  The  fluid 
inlet  temperature  to  this  radiator  is  equal  to  the  outlet 
temperature  of  the  heat  exchanger  cooling  the  electronic 
components.  Likewise,  the  fluid  outlet  temperature  of  the 
radiator  is  equal  to  the  inlet  temperature  of  the  heat 
exchanger.  For  cooling  of  the  electronic  components,  the 
temperature  of  the  secondary  radiator  working  fluid  can  not 
exceed  the  temperature  of  the  components  being  cooled  (unless  a 
refrigerator  is  used).  Thus,  operation  of  the  components  at 
higher  temperatures  allows  reduction  of  the  secondary  radiator 
area  in  two  ways.  First,  less  waste  heat  needs  to  be  rejected 
actively  to  maintain  a  higher  payload  temperature,  and  second, 
the  inlet  and  outlet  temperatures  of  the  secondary  radiator  can 

be  higher. 

The  PCU  and  PCU  shell  are  both  modeled  as  cylinders  with  the 
outer  cylinder  representing  the  shell.  The  PCU  consists  of  the 
electronic  components  required  to  convert  the  platform 
electrical  power  into  a  form  which  the  payload  can  use.  These 
components  typically  are  about  95%  efficient  and  thus  a 
relatively  small  amount  of  waste  heat  will  be  generated  within 
the  PCU. 

The  payload  and  its  shell  are  geometrically  modeled  the  same 
way  as  the  PCU  and  its  shell.  The  payload  consists  of  electronic 
equipment  required  for  the  operation  of  the  platform  such  as 
computers,  radars,  radios,  and  guidance  systems.  Essentially  all 
the  electrical  power  consumed  by  these  components  will  be 
converted  to  waste  heat  (with  the  possible  exception  of  the 
radars) .  Active  cooling  of  these  components  will  therefore  be 
essential. 

Modeling  all  the  electrical  equipment  as  a  single  payload  is 
based  on  the  simplification  that  the  different  electronic 
components  operate  at  the  same  temperature.  In  reality,  each 
component  will  have  different  temperature  limitations  and  may 
even  require  its  own  radiator . 

Surrounding  the  PCU  and  payload  will  be  shells.  These  shells 
serve  two  purposes;  (1)  to  provide  protection  of  the  internal 
components  from  meteoroids,  and  (2)  to  provide  additional 
surface  area  for  radiation  of  waste  heat  to  space.  Thus,  the  PCU 
and  payload  each  have  a  "shell  radiator"  surrounding  them. 
Transfer  of  heat  from  the  internals  to  the  shell  will  occur  by 
radiation  and  via  conduction  through  any  support  structure. 
Also,  heat  pipes  (which  have  very  high  effective  conductivities) 
may  be  connected  between  the  shell  and  the  internals  to  take 
maximum  advantage  of  the  additional  shell  surface  area. 

A  computer  program  has  been  written  to  model  the  heat 
transfer  between  the  major  structures  of  this  conceptual  space 
platform.  Because  the  heat  flux  across  the  radiators  will  vary 
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with  position,  the  program  allows  the  user  to  divide  the 
radiators  into  any  number  of  nodes.  Because  view  factors  for  all 
the  surfaces  of  the  platform  cannot  be  determined  until  the 
sizes  of  the  radiators  are  known,  an  iterative  approach  is  used 
to  solve  for  the  radiator  dimensions.  Also,  the  maximum 
temperature  at  which  the  secondary  radiator  can  operate  depends 
on  the  temperatures  of  the  PCU  and  payload.  Again,  an  iterative 
scheme  is  used  to  determine  this  temperature.  If  a  refrigerator 
is  used  to  cool  the  payload,  the  program  determines  the 
additional  thermal  power  required  to  run  the  refrigerator  along 
with  the  resulting  increase  in  the  main  radiator  size.  Finally, 
heat  transfer  from  the  PCU  or  payload  to  the  PCU  or  payload 
shell  will  occur  both  by  radiation  and  by  conduction.  Coupling 
of  the  two  modes  of  heat  transfer  is  accomplished  by  using  the 
Modified  Regula  Falsi  method  for  finding  roots  of  a  function. 

This  report  contains  a  description  of  the  models  used  in  the 
program  to  determine  the  radiator  sizes  and  the  PCU  and  payload 
temperatures.  Also,  the  results  of  parametric  calculations  are 
presented  to  demonstrate  the  heat  transfer  dependence  on  various 
input  parameters  such  as  the  number  of  nodes,  the  separation 
distances  between  the  various  structures,  the  conduction 
parameters,  and  the  use  of  a  refrigerator. 
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RADIATOR  RADIATOR 


Figure  1.1  Schematic  Diagram  of  Platform  Structures 


2.0  MATHEMATICAL  MODELS 


2.1  Radiosity  Equations 

The  first  step  toward  determination  of  the  radiator  sizes 
and  the  payload  and  PCU  temperatures  is  the  development  of  a 
model  to  calculate  the  radiant  heat  transfer  between  multiple 
surfaces.  For  this  model,  the  assumption  has  been  made  that  all 
the  surfaces  are  gray  and  diffuse.  A  gray  surface  is  one  whose 
radiative  properties  (i.e.,  emissivity)  are  not  a  function  of 
wavelength.  A  diffuse  surface  is  one  that  emits  and  reflects 
radiation  with  equal  intensity  in  all  directions.  The  first 
assumption  is  reasonable  for  this  model  because  most  of  the 
radiation  will  be  emitted  at  relatively  long  wavelengths.  The 
radiation  properties  for  the  surfaces  of  interest  do  not  vary 
much  at  these  long  wavelengths.  The  second  assumption  is  adopted 
because  the  directional  characteristics  of  the  surfaces  simply 
are  not  known.  The  use  of  these  two  assumptions  offers 
considerable  simplification  of  the  radiant  heat  transfer 
equations  and  is  consistent  with  the  level  of  accuracy  required 
keeping  in  mind  the  conceptual  nature  of  the  space  platform. 
(Such  a  platform  has  never  been  built  or  designed.) 

To  model  the  radiant  heat  transfer  between  the  surfaces,  a 
method  known  as  the  radiosity  method  [1]  was  used.  This  method, 
common  in  the  literature,  allows  the  calculation  of  radiant  heat 
transfer  between  any  number  of  surfaces.  The  radiosity  of  a 
surface  is  defined  as  the  rate  at  which  radiation  leaves  a 
surface  by  both  emission  and  reflection.  For  N  surfaces,  a 
linear  set  of  N  equations  with  N  unknowns  is  solved  to  determine 
the  radiosity  for  each  surface.  From  the  radiosities,  the  heat 
flux  or  temperature  can  be  determined  for  the  surfaces. 

Three  types  of  surfaces  are  considered  in  the  radiosity 
method:  (1)  surfaces  for  which  the  temperature  is  known,  (2) 
surfaces  for  which  the  net  heat  flux  is  known,  and  (3) 
"two-sided"  surfaces.  A  two-sided  surface  is  one  in  which  the 
temperatures  of  both  sides  of  the  surface  are  equal,  although 
unknown.  Also,  the  sum  of  the  heat  fluxes  from  both  sides  of  the 
surface  is  known;  however,  the  heat  flux  for  each  side  is 
unknown.  An  example  of  a  two-sided  surface  is  a  structure 
radiating  to  space  (a  radiator)  with  the  sun  shining  on  only  one 
side.  Thus,  the  heat  flux  leaving  each  side  of  the  radiator  is 
unknown  even  though  the  total  heat  produced  by  the  radiator  is 
known.  If  heat  conduction  from  one  side  of  the  radiator  to  the 
other  is  very  large  (i.e.,  low  thermal  resistance),  the  two 
sides  will  have  the  same  temperature;  however,  this  temperature 
is  unknown. 
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The  equations  relating  the  radiosity  for  one  surface  to  the 
radiosity  of  all  the  other  surfaces  are  derived  in  reference  1 
for  surface  types  1  and  2;  i.e.,  surfaces  with  either  known 
temperature  or  net  heat  flux.  These  equations  are  given  below. 


(1) 


(2) 


where:  Ji  =  radiosity  for  surface  i, 

Tj  =  temperature  for  surface  i, 

9^/A.  =  heat  flux  for  surface  i, 

=  emissivity  for  surface  i, 

Pj :  =  view  factor  from  surface  i  to  surface  j, 
oJ  =  Stefan-Boltzmann  constant. 


and 


The  equations  for  surface  type  3,  a  two-sided  surface,  are 
derived  below  starting  with  equations  (1)  and  (2)  along  with  the 
additional  constraints  that: 


=  Tk 


and. 


(3) 


Q  =  qi  +  qk  (4) 

The  subscripts  i  and  k  refer  to  the  two  sides  of  the 
surface.  In  the  matrix  of  radiosity  equations,  each  side  of  a 
two-sided  surface  is  represented  as  a  separate  surface. 
Therefore,  information  indicating  what  surfaces  are  "connected" 
to  form  a  two-sided  surface  must  be  provided  within  the  program. 

First,  equation  (1)  is  written  for  both  sides  of  the 
two-sided  surface,  i  and  k.  These  two  equations  can  then  be 
combined,  using  the  first  constraint  given  above  (equation  3) , 
to  eliminate  the  unknown  temperatures.  Thus, 


J. 


l 


(5) 
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Now,  equation  (2)  is  written  for  both  sides  of  the  surface 
and  the  second  constraint  (equation  4)  is  used  to  eliminate  the 
unknown  heat  fluxes.  Also  note  that  the  areas  of  both  sides  are 
equal  and  recall  that  Q  is  the  total  heat  produced  by  the 
two-sided  surface.  Thus, 


8- 


(6) 


Equations  (5)  and  (6)  are  now  rewritten  such  that  the 
radiosity  for  the  ith  surface  is  isolated  on  the  left  side  of 
the  equal  sign  and  all  other  terms  are  placed  on  the  right. 
(This  is  done  to  place  the  equations  in  a  form  suitable  for  use 
with  a  Gauss-Siedel  iterative  scheme  and  to  simplify  the 
derivation.)  The  new  equations  are: 


These  summations  are  over  all  surfaces  except  the  ith 
surface  because  this  term  has  been  factored  out.  It  is  necessary 
to  satisfy  both  the  temperature  and  the  heat  flux  constraints 
simultaneously  for  both  sides  of  the  two-sided  surface. 
Therefore,  equation  (6')  is  written  for  the  kth  surface  and 
then  substituted  into  equation  (5')  to  yield  a  single  equation 
for  the  radiosity  of  the  ith  surface.  After  simplification  and 
rearrangement : 


Thus,  equations  (1) ,  (2) ,  and  (7)  can  be  used  to  determine 
the  radiosities  for  surfaces  within  an  enclosure  containing 
surface  types  1,  2,  and  3,  respectively.  As  already  mentioned,  a 
Gauss-Siedel  iterative  technique  was  chosen  to  solve  these 
equations.  In  this  technique,  the  initial  radiosities  for  all 
surfaces  are  guessed.  Using  these  initial  values,  new  values  are 
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calculated  for  each  surface  with  the  updated  value  used  for  each 
successive  radiosity  calculation.  With  this  technique,  any 
desired  accuracy  (which  should  be  consistent  with  the  accuracy 
of  the  radiation  properties)  can  be  chosen  for  the  desired 
convergence  criteria. 

These  equations  were  programmed  for  computer  application  in 
a  subroutine  called  RADHT.  Because  of  the  generality  with  which 
the  equations  are  written,  it  is  possible  to  use  this  subroutine 
to  solve  for  the  radiation  heat  transfer  within  any  enclosure 
simply  by  changing  the  input  to  the  subroutine.  To  lend  some 
confidence  in  the  use  of  this  subroutine,  several  assessment 
calculations  were  performed.  The  problems  for  these  calculations 
were  taken  from  various  textbooks  and  are  presented  in  Appendix 
A.  Very  good  agreement  with  the  textbook  solutions  was  achieved 
for  all  cases. 

2.2  Radiator  Sizing  Equations 

The  platform  radiators  have  been  assumed  to  be  contact  heat 
exchangers  as  discussed  in  reference  2.  In  a  contact  heat 
exchanger,  a  working  fluid  flowing  in  conduits  transfers  waste 
heat  to  heat  pipes  which  are  "plugged  into"  the  conduit.  The 
heat  pipes  then  radiate  the  waste  heat  to  space.  For  the  disk 
radiators  in  this  study,  the  conduits  are  like  the  spokes  of  a 
wheel  and  the  heat  pipes  are  connected  circumferentially  to  form 
the  disk.  In  a  contact  radiator,  the  contact  resistance  between 
the  conduit  and  the  heat  pipes  can  be  made  very  small  such  that 
the  temperature  drop  across  the  contact  can  be  considered 
insignificant.  Also,  the  temperature  drop  from  the  working  fluid 
to  the  conduit  wall  is  typically  less  than  5  K  [2].  Thus,  very 
efficient  use  of  the  heat  transfer  area  can  be  achieved.  For  the 
purposes  of  determining  radiator  size,  the  temperature  drop  from 
the  fluid  to  the  heat  pipes  was  assumed  to  be  zero.  Because  this 
temperature  drop  is  really  nonzero,  the  actual  radiator  size 
would  be  somewhat  larger. 

The  input  to  subroutine  RADHT,  discussed  in  the  previous 
section,  includes  the  view  factors  between  all  the  different 
surfaces  of  the  enclosure.  The  problem  arises  that  the  view 
factors  can  not  be  determined  until  the  dimensions  of  the 
radiators  are  known.  To  further  complicate  this  problem  is  the 
fact  that  the  temperature  distribution  from  the  inside  edge  to 
the  outside  edge  of  each  disk  radiator  is  also  not  known.  This 
section  includes  a  discussion  of  an  iterative  approach  to  solve 
this  problem. 

First,  consider  a  radiator  in  the  shape  of  a  disk.  Assuming 
T  =  T-^  =  T2/  an  energy  balance  on  a  differential  area 
element,  dA,  of  the  radiator  yields: 
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where: 


—  surface  emissivity  of  side  1, 
c2  =  surface  emissivity  of  side  2, 

=  surface  temperature  of  side  1, 

T2  =  surface  temperature  of  side  2, 

Tg,  =  effective  background  temperature,  side  1, 

Tg2  =  effective  background  temperature,  side  2, 

©  =  stefan-Boltzmann  constant, 

m  =  working  fluid  mass  flow  rate,  and 
Cp  =  working  fluid  specific  heat. 

The  boundary  conditions  for  this  differential  equation  are  the 
radiator  inlet  and  outlet  temperatures  which  are  program  input 
values.  The  required  value  for  the  quantity  mCp  can  thus  be 
determined  as  the  heat  to  be  rejected  divided  by  the  inlet  to 
outlet  temperature  difference. 

Therefore,  for  steady-state  operation,  the  energy  deposited 
in  the  differential  area  element  by  the  radiator  working  fluid 
must  equal  the  energy  radiated  from  both  sides  of  the  radiator. 
The  solution  of  -.this  equation  yields  the  temperature 
distribution  across  the  radiator.  The  key  to  solving  this 
equation  is  in  finding  an  effective  background  temperature, 
T  ,  that  accounts  for  the  presence  of  all  the  structures  of 
the  platform. 

This  equation  can  be  simplified  by  making  some  algebraic 
substitutions.  First,  for  the  disk: 

dA  =  2ir  rdr 

where  r  is  the  radius. 


Also, 

c  =  +  c2)/2.0 

C  «  -±c  /4ttco 
P 


<C1  TS1 


E2Ts2 


>/(t 


This  last  substitution  allows  the  replacement  of  two  variables 
by  a  single  variable.  This  implies  that  it  is  possible  to 
determine  a  single  effective  background  temperature  for  both 
sides  of  the  element;  i.e.,  it  is  not  necessary  to  determine  the 
effective  background  temperature  for  each  side  separately.  This 
substitution  considerably  simplifies  the  solution  because  fewer 
variables  are  required  and  less  computation  per  iteration  is 
needed.  Also,  an  analytical  solution  to  the  integral  can  be 
found.  Without  this  substitution,  the  integration  would  have  to 
be  performed  numerically. 
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These  substitutions  result  in  the  following  equation: 


The  solution  to  this  differential  equation  can  be  determined 
by  the  integration  of  each  side  of  the  equation  over  the 
appropriate  intervals.  Thus, 


where:  R-r  =  the  disk  inner  radius, 

Rq  =  the  disk  outer  radius, 

Tj  =  the  temperature  at  the  inner  radius,  and 
Tq  =  the  temperature  at  the  outer  radius. 

Performing  the  integration  yields: 


If  the  disk  is  subdivided  into  N  nodes,  equation  (11)  can  be 
applied  to  each  node  where  the  subscripts  i  and  i+1  refer  to  the 
inner  and  outer  boundary  of  the  node,  respectively. 

At  this  point,  two  different  solution  approaches  could  be 
used.  The  first  approach  is  to  use  equation  (11)  to  determine 
the  outer  radius  of  the  radiator  for  a  given  value  of  Ts.  (The 
inner  radius,  inside  temperature,  and  outside  temperature  are 
known.)  Now  the  radiator  can  be  divided  into  N  nodes  of  any 
arbitrary  width.  Equation  (11)  can  then  be  used  to  solve  for  the 
temperatures  at  the  boundaries  of  the  nodes.  These  temperatures 
must  be  solved  iteratively  due  to  the  transcendental  nature  of 
the  equation.  A  simpler  approach  (the  approach  used  in  the 
program)  is  to  assume  that  the  node  widths  will  be  such  that  the 
temperature  drop  across  all  nodes  is  equal.  Thus, 

Ti+1  =  Ti  +  (Tj  -  T0)/N 

Now,  equation  (11)  can  be  used  to  solve  for  the  node  radii 
directly.  This  approach  has  the  additional  benefit  of 
automatically  placing  more  nodes  in  the  region  of  the  largest 
temperature  gradient. 
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With  the  node  radii  and  end  temperatures  thus  computed,  it 
is  necessary  to  determine  an  "average"  temperature  for  each 
node.  This  average  is  not  the  algebraic  average  of  the  node  end 
temperatures  because  of  the  fourth  power  temperature  dependance 
for  radiant  heat  transfer.  To  determine  the  average  node 
temperature,  consider  an  energy  balance  on  the  node* 


(Ti  -  Ti.i) 


2acTT  ri+l 


-  '.KV  -  *.!) 


where  T^  is  the  average  temperature  that  satisfies  this 
equation.  Rearranging  this  equation  yields: 


(12') 


Now  equation  (10)  can  be  used  to  determine  the  radius  at  which 
the  average  temperature  occurs.  To  do  this,  perform  the 
integration  with  r^  as  the  upper  limit  of  integration  for  r 
and  T.^  as  the  upper  limit  of  integration  for  T. 

Thus,  equations  (11)  and  (12')  can  be  used  to  determine  the 
dimensions  and  temperature  distribution  of  the  radiator  for  a 
given  effective  background  temperature,  Ts  i,  for  each  node. 
The  task  now  is  to  determine  the  Ts  values  for  both  radiators. 
The  following  iterative  procedure  is  used  to  accomplish  this 
task. 

STEP  1:  Use  the  effective  temperature  of  space  as  a  first  guess 
to  Ts  for  each  node. 

STEP  2:  Use  equations  (11)  and  (12')  to  determine  the  node  radii 
and  temperatures. 

STEP  3:  Now  that  values  for  the  radiator  dimensions  exist,  view 
factors  for  all  the  surfaces  of  the  platform  can  be  calculated. 
Also,  values  for  the  average  node  temperatures  exist.  Therefore, 
the  subroutine  RADHT  can  be  used  to  determine  the  heat  fluxes 
and  temperatures  for  all  surfaces. 

STEP  4:  Now  calculate  the  power  radiated  away  by  both  sides  of 
each  node  using  the  RADHT  calculated  values  of  heat  flux.  Thus, 

Pi  =  qi  +  qjc  (13> 

where  the  q's  are  the  product  of  the  node  heat  flux  and  area. 
The  i  and  k  subscripts  indicate  the  two  sides  of  the  node. 

Now,  for  each  node, 

PA  =  2Ai al  (Ti4  -  TSfi4)  (14) 
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Rearranging  yields  an  expression  for  updating  the  values  of 

TS' 

Ts,i'  =  (Ti4  -  Pi/A ±2  oe  )1/4  (14') 

STEP  5:  Now,  the  updated  values  of  Ts  can  be  compared  to  the 
previous  values.  Steps  1  through  4  can  then  be  repeated  until 
these  values  converge  or  when  the  Pi  values  converge  to  the 
desired  degree  of  accuracy. 

Use  of  the  Gauss-Siedel  iterative  technique  for  solving  the 
radiosity  equations  fits  in  well  with  this  procedure  because  the 
radiosities  from  the  previous  iteration  can  be  used  for  the 
initial  guess  for  the  current  iteration. 

Subroutine  RADSIZE  contains  the  logic  to  perform  the 
radiator-sizing  calculations.  The  subroutine  has  been  written  to 
allow  the  radiator  inlet  temperature  to  be  at  either  the  inside 
or  outside  radius  of  the  radiator.  As  will  be  discussed  later, 
it  may  be  advantageous  (with  respect  to  minimizing  required 
radiator  area)  to  deliver  the  hot  working  fluid  to  the  outside 
of  the  radiator  first,  via  insulated  conduits,  and  to  then 
collect  the  cooler  fluid  at  the  inner  radius. 

If  the  radiator  is  to  reject  heat  at  a  constant  temperature 
(as  would  be  the  case  for  a  Rankine  cycle  main  power  supply) , 
the  radiator  inlet  and  outlet  temperatures  are  equal  and  thus 
the  temperature  distribution  across  the  radiator  is  known.  For 
this  case,  the  preceding  procedure  is  simplified.  Now,  equation 
(14)  can  be  used  to  solve  for  the  node  areas  directly  because 
the  node  temperatures  are  known.  However,  the  effective 
background  node  temperatures  remain  as  unknowns  and  the 
iterative  approach  is  still  required. 

2.3  View  Factors 

It  is  necessary  to  calculate  view  factors  between  all 

surfaces  of  the  platform.  This  is  done  in  subroutines  VIEW1  and 
VIEW2.  Subroutine  VIEW1  determines  the  view  factors  between  the 
payload  and  PCU  and  their  shells.  The  view  factors  between  these 
structures  need  to  be  calculated  only  once  because  the 
dimensions  for  these  structures  are  fixed.  Subroutine  VIEW2 

determines  the  view  factors  between  all  the  nodes  of  the  two 

radiators  and  all  other  surfaces.  Because  the  dimensions  of  the 
radiators  are  determined  iteratively,  it  is  necessary  to 

recalculate  these  view  factors  for  each  iteration. 

Both  view  factor  subroutines  make  use  of  only  two  basic  view 
factor  expressions.  All  other  view  factors  are  determined  by 
view  factor  algebra  or  by  using  the  fact  that  the  sum  of  the 
view  factors  from  one  surface  to  all  the  others  must  equal  1.0.  ' 
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The  two  basic  view  factor  expressions  are  for  a  disk-to-disk  and 
a  cylinder-to-concentric-cylinder.  The  expressions  for  the  two 
view  factors  are  taken  from  reference  3 .  A  considerable  amount 
of  view  factor  algebra  is  required  to  arrive  at  expressions  for 
view  factors  between  cylinders  and  rings  and  between  rings  and 
other  rings.  Figure  2.3.1  provides  schematic  diagrams  of  the 
different  geometries  for  which  view  factors  are  calculated. 
Subroutines  VIEW1  and  VIEW2  make  use  of  other  subroutines  and 
functions  to  calculate  all  the  required  view  factors.  The 
following  is  a  list  of  these  subroutines  and  functions. 

FUNCTION  DTOD  -  determines  the  view  factor  between  a  disk  of 
radius  r-j^  to  another  disk  of  radius  r2  separated  by  distance 
h.  (See  reference  3,  page  826,  #21.) 

FUNCTION  CTOD  -  determines  the  view  factor  between  a  cylinder  of 
radius  r-j^  and  an  adjacent  and  perpendicular  disk  of  radius 
r^.  Some  view  factor  algebra  is  included  starting  with  the 
view  factor  between  two  concentric  cylinders.  (See  reference  3, 
page  828,  #28.) 

SUBROUTINE  FRINGS  -  uses  function  DTOD  and  view  factor  algebra 
to  determine  ring-to-ring  view  factors.  The  rings  are  separated 
by  distance  h. 

SUBROUTINE  CTORING  -  uses  function  CTOD  and  view  factor  algebra 
to  determine  cylinder-to-ring  view  factors.  The  cylinder  and 
ring  are  separated  by  distance  h. 

Because  the  radiators  can  be  divided  into  any  number  of 
nodes,  the  view  factor  subroutines  must  be  written  with  the 
number  of  nodes  as  a  parameter.  The  details  of  the  bookkeeping 
and  the  view  factor  algebra  used  to  arrive  at  all  the  view 
factors  will  not  be  explained  in  any  further  detail.  However, 
Appendix  B,  which  contains  a  computer  listing  of  the  entire 
program,  can  be  examined  for  the  appropriate  FORTRAN 
expressions. 
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Figure  2.3.1  View  Factor  Geometries 


2.4  Coupled  Heat  Transfer  Modes 

The  radiator-sizing  equations  discussed  in  the  previous 
sections  can  be  solved  without  considering  conduction  between 
the  payload  (or  PCU)  internals  and  the  shell  because  only  the 
total  heat  flux  leaving  the  shell  is  of  consequence  for  radiator 
sizing.  In  other  words,  it  doesn't  matter  if  the  waste  heat 
generated  in  the  internals  is  transferred  to  the  shell  by 
conduction  or  radiation;  it  is  the  total  heat  transferred  that 
determines  the  shell  temperature  and  heat  flux.  Therefore,  the 
radiator-sizing  calculations  can  be  performed  assuming  that  all 
the  waste  heat  is  transferred  from  the  internals  to  the  shell  by 
radiation  and  that  none  of  the  waste  heat  is  transferred  to  the 
shell  by  conduction.  However,  determining  the  temperature  of  the 
internals  does  depend  on  the  amount  of  waste  heat  removed  by 
conduction  and  radiation. 

Waste  heat  generated  in  the  PCU  and  payload  internals  will 
be  transferred  from  the  internals  to  the  shell  by  both  radiation 
and  conduction.  The  conduction  will  occur  along  the  support 
structure  used  to  support  the  shell.  Also,  heat  pipes  may  be 
connected  between  the  internals  and  shell  to  greatly  improve 
this  heat  transfer  mechanism.  Heat  pipes  have  very  large 
effective  thermal  conductivities  and  essentially  can  transfer 
heat  isothermal ly.  The  coupling  of  the  conduction  and  radiation 
heat  transfer  modes  is  nonlinear  because  conduction  is 
proportional  to  the  first  power  of  temperature;  whereas, 
radiation  is  proportional  to  the  fourth  power.  To  accomplish 
this  coupling,  use  was  made  of  the  Modified  Regula  Falsi  [4] 
method  for  finding  roots  of  a  function.  This  coupling  approach 
is  discussed  in  this  section. 

An  energy  balance  on  the  internals  can  be  written  as: 

q  -  qv  -  <3d  "  Sr  =  °*°  (15) 

where:  q  =  total  waste  heat  generated, 

qv  =  waste  heat  removed  by  convection, 
qd  =  waste  heat  removed  by  conduction,  and 
qr  =  waste  heat  removed  by  radiation. 

The  expression  for  the  convection  term  depends  on  the  type 
of  heat  exchanger  used  to  cool  the  internals;  i.e.,  the  number 
of  fluid  passes,  the  size  and  type  of  fins,  etc.  Given  the 
conceptual  nature  of  the  platform,  convection  heat  transfer  is 
not  modeled  in  detail.  Instead,  the  convection  term  is  treated 
parametrically  such  that  qv  is  assumed  to  be  a  known  quantity. 
(The  convection  term  will  be  discussed  in  more  detail  later  in 
this  section.)  Thus,  two  terms  in  equation  (15)  are  unknown, 
namely,  qd  and  qr.  Equation  (15)  can  be  rewritten  in  the 
form  of  a  function  as: 

f(qr)  =  q  -  qv  -  Sd  "  3r  (15/) 
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Finding  the  values  of  qr  and  qd  that  make  this  function 
equal  to  zero  is  desired.  This  is  equivalent  to  finding  the  root 
of  an  equation  for  which  the  Modified  Regula  Falsi  method  is 
applicable. 

To  evaluate  the  radiation  term,  the  radiosity  equations  of 
Section  2.1  can  be  used  for  the  internals  and  shell  surfaces. 
Solution  of  the  equations  yields  the  shell  and  internal 
temperatures  from  which  the  conduction  term  can  be  calculated 
using: 

qd  =  (kA/x)  AT 

where:  k  =  effective  thermal  conductivity, 

A  =  heat  transfer  area  for  conduction, 

x  =  conduction  path  length,  and 

AT  =  temperature  drop  from  internals  to  shell. 

Thus,  qd  is  a  function  of  the  unknown  internals 
temperature  which  is  a  function  of  qr.  Choosing  qr  =0.0 
(i.e.,  all  the  waste  heat  is  removed  from  the  internals  by 

conduction)  and  qr  =  q  -  qv  (i.e.,  all  the  waste  heat  is 

removed  by  radiation)  provides  an  interval  containing  a  root 
that  satisfies  equation  (15').  Starting  with  this  interval,  the 
Modified  Regula  Falsi  method  quickly  finds  the  root  and  thus  the 
internals  temperature.  This  is  equivalent  to  solving  the 
radiation  and  conduction  equations  simultaneously.  However,  use 
of  the  Modified  Regula  Falsi  method  allows  one  to  keep  the 
expressions  for  conduction  and  radiation  separate.  A  graphical 
depiction  of  this  process  is  shown  in  Figure  2.4.1.  The  program 
was  written  such  that  either  the  kA/x  value  or  the  at  value  can 

be  input.  If  at  is  input,  the  required  kA/x  value  to  achieve 

this  AT  is  calculated. 

Returning  to  the  convection  term,  recall  that  qv  was 
treated  as  a  known  quantity.  To  obtain  an  order  of  magnitude 
estimate  of  the  hA  value  required  to  remove  the  specified  amount 
of  waste  heat,  Newton's  law  of  cooling  can  be  used.  The  waste 
heat  removed  by  convection  is  thus  expressed  as: 


qv  =  hA(T  -  Tc)  (16) 

where:  h  =  effective  heat  transfer  coefficient, 

A  =  total  heat  transfer  area  available, 

T  =  temperature  of  the  internals,  and 
Tc  =  average  coolant  temperature. 

An  estimate  of  the  value  of  Tp  is  made  using  the  algebraic 
average  of  the  heat  exchanger  inlet  and  outlet  temperatures. 
Note  that  for  heat  transfer  from  the  internals  to  the  fluid,  the 
fluid  temperature  must  be  less  than  the  internals  temperature. 
Also,  recall  that  the  outlet  temperature  of  the  heat  exchanger 
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is  equal  to  the  inlet  temperature  of  the  secondary  radiator.  The 
inlet  temperature  of  the  secondary  radiator  must  be  known  before 
the  radiator  can  be  sized.  Therefore,  a  guess  is  made  for  this 
temperature  and  the  radiator-sizing  calculations  and 
conduction-convection  coupling  calculations  are  performed.  An 
updated  value  is  then  calculated  based  on  the  calculated 
temperature  of  the  internals.  Thus, 


Ti  = 
T0  = 


T  “  A^min 

Tj  +  AT 


where:  Tj  =  the  secondary  radiator  inlet  temperature 

(equivalent  to  the  heat  exchanger  outlet  temperature) , 
T0  =  the  secondary  radiator  outlet  temperature 
(equivalent  to  the  heat  exchanger  inlet  temperature) , 
AT  =  the  desired  radiator  temperature  difference, 

T  =  the  internals  temperature,  and 

Atmin  =  the  minimum  allowed  temperature  difference 
lk  1  between  the  internals  and  cooling  fluid. 


Figure  2.4.2  graphically  demonstrates  the  relationship  between 
the  various  temperatures.  The  updated  values  of  Tj  and  Tq 
are  then  used  to  repeat  the  radiator  sizing  and  internal 
temperature  calculations. 
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.4.1  The  Root  of  the  Conduction-Radiation  Equation 


HEAT  secondary 

EXCHANGER  RADIATOR 


Figure  2.4.2 


Radiator  and  Internal  Temperature  Relationship 


2.5  Refrigerator  Model 

The  computer  program  has  been  written  so  that  the  use  of  a 
refrigerator  to  cool  the  payload  internals  can  be  investigated. 
The  possible  benefit  of  a  refrigerator  is  that  it  allows  the 
rejection  of  payload  waste  heat  at  a  temperature  higher  than^the 
payload,  thereby  decreasing  the  size  of  the  required  secondary 
radiator.  However,  because  additional  power  is  required  to  run 
the  refrigerator,  more  waste  heat  associated  with  the  mam  power 
supply  will  be  generated  necessitating  a  larger  main  radiator. 

The  refrigerator  coefficient  of  performance  (COP)  is 
calculated  as  some  specified  fraction  of  the  Carnot  COP.  The 
Carnot  COP  (COPc)  is  calculated  as: 


COPc  =  T-l  /  (Th  -  Tx) 


(17) 


where:  T^  =  the  lowest  temperature  of  the  working  fluid, 

and  ,  .  . 

Th  =  the  highest  temperature  of  the  working  fluid. 

The  secondary  radiator  is  assumed  to  operate  at  Th.  The 
workinq  fluid  is  assumed  to  operate  at  the  payload  temperature 
minus  some  specified  payload-to-working  fluid  temperature 
difference.  Because  the  payload  temperature  is  unknown,  an 
iterative  procedure  is  employed  to  determine  the  value  of  T-^ 
such  that  T-i  equals  the  calculated  payload  temperature  minus 
the  specif ied  payload-to-working  fluid  temperature  difference. 

The  amount  of  electric  power  needed  to  power  the 
refrigerator,  W,  is  determined  as  Q  divided  by  COP  where  Q  is 
the  amount  of  waste  heat  to  be  removed  by  the  refrigerator.  (The 
waste  heat  associated  with  the  refrigerator  is  not  included  in 
the  payload  waste  heat.)  The  amount  of  additional  thermal  power 
required  is  found  as  W  divided  by  the  thermal-to-electric  power 
conversion  efficiency.  An  increase  in  the  power  supply  increases 
the  amount  of  waste  heat  to  be  rejected  by  the  mam  radiator. 
Thus,  the  required  size  of  the  main  radiator  must  be  updated 

each  iteration. 
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3.0  COMPUTER  PROGRAM  APPLICATION 


3.1  Example  Input  and  Output 

To  demonstrate  the  use  and  capabilities  of  the  computer 
program,  several  parametric  calculations  were  performed.  This 
section  includes  a  description  of  these  calculations.  An 
annotated  example  input  file  is  given  in  Table  3.1.1.  The 
annotations  provide  a  concise  description  of  the  required  input 
variables.  This  example  input  provides  a  "base  case"  on  which 
all  subsequent  parametric  calculations  are  based.  The  resulting 
program  output  for  this  example  problem  is  provided  in  Appendix 
C.  Only  three  nodes  for  the  main  radiator  and  four  nodes  for  the 
secondary  radiator  were  used  to  keep  the  output  relatively 
short. 

Table  3.1.1  Sample  Program  Input 


15 

0.01 

20 

0.005 

3 

4 

0.90 

0.90 

0.90 

0.90 

0.9. 0.9. 0.9 
0.9. 0.9. 0.9 
0.9. 0.9. 0.9 
0.9. 0.9. 0.9 
0.80 
0.80 
300.0 
480.0 
965.0 
470.0 
400.0 
6.0 
7.0 
1.0 
30.0 
1.0 
6.0 
12.0 
10.0 
15.0 
6.0 
10.0 
6.0 
12.0 
40 . 0E6 
0.25 
0.0 
C .  05 
0.0 
0.9 
0.9 
-5.0 
-2.0 
-2.0 


MAXIMUM  ITERATIONS  FOR  RADIATOR  SIZING  (-  FOR  DEBUG) 

RADIATOR-SIZING  RELATIVE  CONVERGENCE  CRITERIA 

MAXIMUM  ITERATIONS  FOR  GAUSS-SIEDEL  SOLUTION 

GAUSS-SIEDEL  RELATIVE  CONVERGENCE  CRITERIA 

NUMBER  OF  MAIN  RADIATOR  NODES 

NUMBER  OF  SECONDARY  RADIATOR  NODES 

EMISSIVITY  OF  LEFT  SIDE  OF  MAIN  RAITAT0R 

EMISSIVITY  OF  RIGHT  SIDE  OF  MAIN  RADIATOR 

EMISSIVITY  OF  LEFT  SIDE  OF  SECONDARY  RADIATOR 

EMISSIVITY  OF  RIGHT  SIDE  OF  SECONDARY  RADIATOR 

EMISSIVITY  OF  INSIDE  OF  PCU  SHELL  (LEFT . CENTER . RIGHT) 

EMISSIVITY  OF  OUTSIDE  OF  PCU  SHELL  (LEFT . CENTER . RIGHT ) 

EMISSIVITY  OF  INSIDE  OF  PAYLOAD  SHELL  (LEFI . CENTER. . RIGHT ) 

EMISSIVITY  OF  OUTSIDE  OF  PAYLOAD  SHELL  (LEFT . CENTER . RIGHT ) 

EMISSIVITY  OF  PCU 

EMISSIVITY  OF  PAYLOAD- 

EFFECTIVE  TEMPERATURE  OF  SPACE  (X) 

TEMPERATURE  AT  INNER  RADIUS  OF  MAIN  RADIATOR  ex'; 

TEMPERATURE  AT  OUTER  RADIUS  OF  MAIN  RADIATOR  (X) 

TEMPERATURE  AT  INNER  RADIUS  OF  SECONDARY  RADIATOR  (X) 

TEMPERATURE  AT  OUTER  RADIUS  OF  SECONDARY  RADIATOR  (X) 

INNER  RADIUS  OF  MAIN  RADIATOR  (M) 

INNER  RADIUS  OF  SECONDARY  RADIATOR  (M) 

SEPARATION  DISTANCE  BETWEEN  PCU  AND  MAIN  RADIATOR  CM) 
SEPARATION  DISTANCE  BETWEEN  MAIN  AND  SEC  RADIATOR  CM) 
SEPARATION  DISTANCE  BETWEEN  SEC  RADIATOR  AND  PAYLOAD  (M) 
DIAMETER  OF  PCU  SHELL  (M) 

LENGTH  OF  PCU  SHELL  (M) 

DIAMETER  OF  PAYLOAD  SHELL  (M) 

LENGTH  OF  PAYLOAD  SHELL  (M) 

DIAMETER  OF  PCU  (M) 

LENGTH  OF  PCU  (M) 

DIAMETER  OF  PAYLOAD  (K) 

LENGTH  OF  PAYLAOD  (M) 

TOTAL  FLATFORM  THERMAL  POWER  (V) 

TKERMAL-TO-ELECTRIC  CONVERSION  EFFICIENCY 

FRACTION  OF  CARNOT  REFRIGERATOR  COP  (IF  0.0.  NO  REFF IGERATO: 
FRACTION  OF  ELECTRIC  POWER  CONVERTED  TO  PCU  WASTE  HEAT 
FRACTION  OF  PCU  WASTE  HEAT  REMOVED  ACTIVELY  (TO  SEC  RAD) 
FRACTION  OF  ELECTRIC  POWER  CONVERTED  TO  PAYLOAD  WASTE  HEAT 
FRACTION  OF  FAYLOAD  WASTE  HEAT  REMOVED  ACTIVELY  (SEC  RAD) 
MIN  DELTA-T  BETWEEN  PAYLOAD  AND  SEC  RAD  WORKING  FLUID  * 

KA  X  (-)  OR  DESIRED  DELTA-T  (-)  FOR  PCU. SHELL  CONDUCTION 
XA  X  (-)  OR  DESIRED  DELTA-T  (-)  FOR  PAYLOAD  SHELL  CONDUCT! 0) 


*  IF  ■  0,  VARIABLE  NOT  USED  AND  THE  INPUT  VALUES  FOR  SECONDARY  RADIATOR 
INLET  AND  OUTLET  TEMPERATURE  ARE  USED  AS  INPUT  (NO  ITERATION). 

IF  A  REFRIGERATOR  IS  USED,  THIS  IS  THE  TEMPERATURE  DIFFERENCE  BETWEEN 
THE  PAYLOAD  AND  THE  REFRIGERATOR  WORKING  FLUID. 
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The  first  page  of  the  output  provides  information  concerning 
the  nodes  of  the  two  radiators,  such  as  the  node  radius,  area, 
temperature,  and  the  effective  background  temperature  for  each 
node.  This  information  is  calculated  by  subroutine  RADSIZE.  The 
next  page  provides  the  final  output,  such  as  radiator  areas,  and 
PCU  and  payload  temperatures.  The  third  page  provides  surface 
identification  numbers  for  all  the  surfaces  used  in  the  radiant 
heat  transfer  solution  as  computed  in  subroutine  RADHT.  Next  is 
a  list  of  the  view  factors  between  all  the  surfaces.  Even  with 
only  a  few  nodes  for  each  radiator,  there  are  841  view  factors 
for  this  example  problem.  Many  of  these  view  factors  are  zero, 
however,  and  are  not  printed.  The  last  page  of  output  displays 
the  results  of  the  final  call  to  subroutine  RADHT;  this  includes 
the  radiosities  for  each  of  the  29  surfaces. 

3.2  Radiator  Nodes 

The  number  of  nodes  used  for  each  radiator  affects  the  final 
calculated  radiator  area  with  the  greater  number  of  nodes 
providing  the  more  accurate  results.  Table  3.2.1  shows  the 
effect  of  varying  the  number  of  nodes  in  the  secondary  radiator 
(the  number  of  nodes  in  the  main  radiator  is  held  constant)  for 
two  different  values  of  waste  heat  delivered  to  the  secondary 
radiator.  This  demonstrates  that  the  choice  of  the  number  of 
nodes  depends  on  the  particular  problem  being  solved;  thus,  care 
should  be  taken  in  selecting  an  appropriate  number  of  radiator 
nodes.  (The  separation  distance  between  the  two  radiators  was 
also  found  to  have  a  big  influence  on  the  number  of  nodes 
required. ) 

Table  3.2.1  Effect  of  the  Number  of  Radiator  Nodes 


Waste  Heat  =  8.1  MW 


# 

Nodes 

Secondary 

Radiator  Area  (m2) 

3 

6198 

6 

6093 

12 

6078 

20 

6075 

Waste  Heat  =  4 

.5  MW 

# 

Nodes 

Secondary 

Radiator  Area  (m2) 

6 

3905 

12 

3042 

20 

2976 

30 

2966 

Figure 

3.2.1 

demonstrates 

the  convergence  of  the  effective 

background 

temperature  used  in 

determining  radiator  size.  Twelve 

nodes  were  used  in  the  secondary  radiator  for  this  example.  The 
initial  guess  for  the  effective  background  temperature  for  each 
node  was  300  K.  The  values  of  the  twelve  node  temperatures  for 
the  next  three  iterations  (the  curves  labeled  2,  3,  and  4)  show 
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that  convergence  occurs  after  some  oscillations.  Originally,  the 
oscillations  were  even  more  severe.  To  dampen  the  oscillations, 
the  updated  Ts  values  were  modified  according  to: 

Ts  =  w  Ts (current)  +  (1-w)  Ts (previous) 

where  w  was  chosen  to  be  0.8.  This  relaxation  value,  determined 
by  trial  and  error,  was  found  to  work  reasonably  well  for  most 

problems . 

Table  3.2.2  shows  information  related  to  the  iterations 
reguired  to  determine  the  radiator  sizes.  In  this  table,  L  is 
the  radiator-sizing  iteration  number  and  M  is  the  number  of 
iterations  required  for  the  radiosity  equations  (subroutine 
RADHT)  to  converge.  The  total  power  is  the  sum  of  all  the  node 
powers.  For  this  example,  the  desired  main  radiator  total  power 
was  30.0  MW  and  the  desired  secondary  radiator  total  power  was 
8.1  MW.  Thus,  the  radiator  sizing  iterative  scheme  converges 
rapidly  for  this  problem. 


Table  3.2.2  Iteration  Summary 


L 


M  Maximum  Node  Power  Total  Power 
Relative  Error  (main  radiator) 


Total  Power 
(sec  radiator) 


1  9 

2  4 

3  2 

4  1 


0.7261 

0.3110 

0.1309 

0.0074 


28.71  MW 
29.16  MW 
30.00  MW 
30.00  MW 


2.489  MW 
10.450  MW 
7.864  MW 
8.090  MW 


All  of  the  calculations  discussed  in  this  section  are  based 
on  the  use  of  an  effective  temperature  of  space  equal  to  300  K. 
(The  effective  temperature  of  space  should  not  be  confused  with 
the  effective  background  temperature  which  accounts  for  the 
presence  of  all  platform  structures  and  space.)  The  effective 
space  temperature  accounts  for  the  radiation  arriving  from  the 
sun  (directly  or  reflected  off  the  earth)  along  with  the 
radiation  emitted  by  the  earth.  This  value  is  a  function  of 
platform  orientation,  altitude,  orbit  inclination,  geometry,  and 
surface-coating  properties  and  can  be  expected  to  range  anywhere 
from  200  K  to  360  K  for  anticipated  platform  applications.  The 
value  of  300  K  was  chosen  as  a  representative  intermediate 

value. 


Figures  3.2.2  and  3.2.3  show  the  main  and  secondary  radiator 
temperatures,  respectively.  The  main  radiator  has  10  nodes  and 
the  secondary  radiator  has  12  nodes.  Also,  the  radiators  are  30 
m  apart.  The  fluid  inlet  to  the  main  radiator  was  chosen  to  be 
at  the  outer  radius  of  the  disk  while  the  inlet  for  the 
secondary  radiator  is  at  the  inner  radius.  This  arrangement 
resulted  in  the  smallest  total  radiator  area.  If  the  inlet  to 
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the  secondary  radiator  is  reversed,  the  cooler  nodes  close  to 
the  inner  radius  do  not  reject  heat  very  well  because  the 
presence  of  the  hot  main  radiator  and  payload  shell  result  in  a 
very  high  effective  background  temperature  for  those  nodes.  An 
alternative  to  reversing  the  flow  direction  would  be  to  increase 
the  inner  radius  of  the  secondary  radiator,  thereby  effectively 
moving  the  radiator  away  from  the  other  platform  structures. 

The  effective  background  temperature  for  the  main  radiator 
nodes  is  relatively  flat  indicating  that  the  other  surfaces  of 
the  platform  do  not  have  a  significant  impact  on  the  main 
radiators  ability  to  reject  heat.  This  is  because  the  main 
radiator  operates  at  a  high  temperature  compared  to  other 
surfaces.  The  effective  background  temperature  for  the  secondary 
radiator  nodes,  however,  is  large  relative  to  the  radiator 
temperature  because  the  temperatures  of  the  other  surfaces  of 
the  platform  are  relatively  high.  These  two  figures  also  show 
that  the  temperature  distributions  of  the  two  radiators  are 
completely  different.  This  is  due  to  the  different  relative 
influence  of  each  platform  surface  on  each  radiator  because  of 
their  different  spacial  orientations.  Changing  the  flow 
direction,  the  PCU  and  payload  sizes,  and  the  platform  structure 
separation  distances  all  have  significant  impact  on  the  radiator 
temperature  distribution. 
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TEMPERATURE  (K>  X10°  EFFECTIVE  BACKGROUND  TEMPERATURE  (K) 


Figure  3.2.1  Radiator  Sizing  Iterations 


Figure  3.2.2  Main  Radiator  Temperature  Distribution 
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TEMPERATURE  (K) 


<480 


Figure  3.2.3  Secondary  Radiator  Temperature  Distribution 
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3.3  Conduction 

Figure  3.3.1  shows  the  effect  of  payload-to-shell  conduction 
on  the  payload  temperature.  It  is  not  clear  if  sufficient  room 
would  exist  for  enough  heat  pipes  to  achieve  the  kA/x  values 
shown  on  this  curve.  However,  a  set  of  parametric  calculations 
was  performed  to  show  the  possible  benefit  of  using  heat  pipes 
to  remove  payload  waste  heat.  The  payload  temperature  drops  very 
rapidly  as  the  kA/x  value  is  increased.  It  is  clear  that  without 
any  heat  pipes  the  payload  temperature  would  be  much  higher. 
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PAYLOAD  TEMPERATURE  (K) 
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Figure  3.3.1  Payload-to-Shell  Conduction 
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3.4  Separation  Distance 

The  geometric  parameter  that  has  the  most  effect  on  the  size 
of  the  radiators  is  the  separation  distance  between  the  main  and 
secondary  radiators.  This  effect  is  demonstrated  in  Figure  3.4.1 
which  shows  the  secondary  radiator  area  as  a  function  of 
separation  distance  for  two  different  amounts  of  waste  heat 
delivered  to  the  secondary  radiator.  It  is  apparent  that 
increasing  the  separation  distance  greatly  reduces  the  required 
size  of  the  radiator.  The  size  of  the  main  radiator  (not  shown 
on  the  Figure)  varies  between  2300  m2  and  1550  m2  for 
separation  distances  from  10  m  to  300  m.  The  effect  of 
separation  distance  on  the  main  radiator  is  not  as  significant 
because  it  operates  at  a  much  higher  temperature  than  the 
secondary  radiator.  However,  if  the  operating  temperature  of  the 
secondary  radiator  is  increased  or  the  temperature  of  the  main 
radiator  decreased,  the  effect  on  the  main  radiator  size  would 
become  more  significant. 

The  radiator  separation  distance  also  influences  the  PCU  and 
payload  temperatures  as  shown  in  Figure  3.4.2.  The  payload 
temperature  is  affected  more  because  the  hot  main  radiator  has 
the  largest  influence  on  the  payloads  ability  to  radiate  waste 
heat  to  space. 

Another  parameter  that  affects  the  size  of  the  radiators  and 
the  payload  and  PCU  temperatures  is  the  separation  distance 
between  the  shells  (PCU  and  payload)  and  the  radiators.  (To  show 
this  effect,  a  radiator-to-radiator  separation  distance  of  500  m 
was  used;  thus,  the  radiators  are  far  enough  apart  to  preclude 
any  influence  on  each  other.)  Figure  3.4.3  shows  the  effect  of 
varying  these  separation  distances  on  the  radiator  areas.  The 
reason  the  radiator  areas  increase  for  separation  distances  from 
about  2  m  to  10  m  is  that  more  of  the  shell  views  the  adjacent 
radiator  in  this  range  of  separation  distances.  At  zero 
separation  distance,  only  the  curved  surface  of  the  shell  views 
the  adjacent  radiator.  As  the  separation  distance  is  increased, 
one  end  of  the  shell  also  begins  to  view  the  radiator,  resulting 
in  an  overall  increase  in  the  shell-to-radiator  view.  Further 
increases  in  the  separation  distance  result  in  an  overall 
decrease  of  this  view  such  that  the  required  radiator  area 
decreases.  A  similar  effect  on  payload  and  PCU  temperatures  is 
shown  in  Figure  3.4.4. 
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Figure  3.4.4  Temperature  Versus  Shell”to—Radiator  Distance 
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3.5  Convection 

One  of  the  output  variables  of  the  computer  program  is  an 
estimate  of  the  value  of  hA  required  for  convective  cooling  of 
tte  Svload.  For  a  given  working  fluid  and  flow  rate,  an 
estimate  of  h  can  be  made;  using  the  estimated  hA  and  h  values, 
the  heat  transfer  area  can  be  determined.  This  crude  estimat 
can  at  least  be  used  to  provide  some  idea  of  the  area  require 
by  the  heat  exchanger  to  actively  cool  the  payload.  Figure J^ad 
•shows  the  effect  of  different  hA  values  on  the  payloa 
temperature.  Increasing  the  value  of  hA  has  a  diminishing 
influence  on  the  temperature  because  the  temperature  difference 
between  the  working  fluid  and  the  payload  decreases  as  more  and 
more  heat  is  removed  from  the  payload. 
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Figure  3.5.1  Effect  of  Convection  on  Payload  Temperature 


3.6  Refrigerators 

The  effect  of  a  refrigerator  was  determined  by  repeating  the 
base  case,  specifying  a  refrigerator  COP  of  70%  of  the  Carnot 
COP.  Also,  the  secondary  radiator  was  assumed  to  operate  at  700 
K  and  the  temperature  difference  between  the  payload  and  the 
refrigerator  working  fluid  was  taken  as  40  K.  (In  the  base  case 
without  a  refrigerator,  the  log  mean  temperature  difference 
between  the  payload  and  the  payload  heat  exchanger  working  fluid 
was  also  approximately  40  K.)  The  results  indicate  that  a  6  MWe 
refrigerator  is  needed  to  cool  the  payload,  requiring  an 
additional  24  MW  of  thermal  power.  Compared  to  the  base  case, 
the  size  of  the  secondary  radiator  decreased  from  6087  irr  to 
395  m^  while  the  main  radiator  increased  in  size  from  1906 
m^  to  2967  m  . 

3.7  Payload  Temperature 

The  principal  motivation  for  creating  this  radiant  heat 
transfer  computer  program  was  to  estimate  the  reduction  in 
secondary  radiator  area  that  could  be  achieved  by  allowing 
higher  payload  and  PCU  operating  temperatures.  Electronic 
components  capable  of  high  temperature  operation  require  less 
active  cooling.  Also,  higher  component  temperatures  allow 
rejection  of  waste  heat  at  higher  temperature  by  the  secondary 
radiator.  The  relationship  between  the  secondary  radiator  area 
and  payload  temperature  is  provided  in  Figure  3.7.1  for  the  base 
case  input  parameters.  This  curve  shows  that  increasing  the 
permissible  payload  temperature  from  450  K  to  600  K  results  in 
over  a  factor  of  two  reduction  in  the  radiator  area. 
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Figure  3.7.1  Effect  of  Payload  Temperature  on  Radiator  Area 
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4.0  SUMMARY  AND  CONCLUSIONS 

A  conceptual  space  platform  consists  of  a  payload,  a  PCU, 
and  two  radiators:  the  main  radiator  and  a  secondary  radiator.  A 
computer  program  was  written  to  determine  the  required  size  of 
the  two  radiators  and  the  temperatures  of  the  PCU  and  payload 
for  a  given  platform  power  level.  The  program  includes  a  model 
for  the  radiant  heat  transfer  between  the  various  platform 
surfaces.  An  iterative  algorithm  is  employed  in  conjunction  with 
this  model  to  determine  the  size  of  the  radiators.  The  program 
user  can  subdivide  the  two  radiators  into  any  number  of  nodes  to 
increase  the  accuracy  of  the  radiant  heat  transfer  solution.  The 
use  of  more  nodes  also  allows  better  prediction  of  the  nonlinear 
temperature  drop  that  occurs  across  the  radiators  as  the  working 
fluid  deposits  the  platform's  waste  heat  in  the  radiator.  View 
factor  expressions  are  automatically  calculated  for  different 
choices  of  the  number  of  nodes.  The  user  can  also  select 
different  separation  distances  between  the  various  platform 
structures.  A  model  is  included  to  couple  the  radiant  and 
conduction  heat  transfer  that  occurs  between  the  payload  and  its 
shell  and  between  the  PCU  and  its  shell.  Also,  the  program 
allows  the  use  of  a  refrigerator  to  cool  the  payload.  If  a 
refrigerator  is  used,  the  program  determines  the  amount  of 
additional  thermal  power  needed  to  run  the  refrigerator. 

This  computer  program  was  used  to  perform  a  variety  of 
parametric  calculations.  The  results  indicate  that  the  secondary 
radiator  size  is  strongly  dependent  on  its  proximity  to  the  main 
radiator  and  that  the  radiant  interchange  that  occurs  between 
the  various  platform  structures  has  a  large  effect  on  the 
structure  temperatures.  The  size  of  the  secondary  radiator  also 
depends  on  the  amount  of  waste  heat  that  it  must  radiate  to 
space  and  the  temperature  at  which  it  operates.  Use  of  a 
refrigerator  to  cool  the  payload  significantly  reduces  the  size 
of  the  secondary  radiator  but  at  the  expense  of  increased  main 
radiator  and  power  supply  size.  Results  of  this  computer  program 
indicate  that  with  or  without  a  refrigerator,  the  use  of  high 
temperature  electrical  components  in  the  PCU  and  payload  can 
significantly  reduce  the  required  size  of  the  secondary  radiator 
by  allowing  a  reduction  in  the  amount  of  waste  heat  to  be 
removed  and  by  allowing  an  increase  in  the  inlet  and  outlet 
temperatures  of  the  radiator. 

The  calculations  in  this  report  were  included  to  demonstrate 
the  intended  use  of  the  program.  In  order  to  perform  these 
calculations  for  a  conceptual  space  platform  it  was  necessary  to 
assume  values  for  many  of  the  program  input  variables.  It  should 
be  noted  that  there  are  a  lot  of  variables  to  consider  and  that 
the  results  can  only  indicate  trends  in  the  various  functional 
relationships.  Assuming  representative  variable  values  and 
performing  parametric  calculations  is  the  best  one  can  do  until 
a  detailed  platform  design  has  been  completed. 
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APPENDIX  A  -  ASSESSMENT  PROBLEMS 


SAMPLE  PROBLEM  (1)  -  CYLINDRICAL  CAVITY  AT  CONSTANT 
TEMPERATURE  REJECTING  HEAT  TO  BLACKBODY  SINK  AT  293.3  K 


RADIATION  ANALYSIS 

SURFACE 

EMISSIVITY 

TEMPERATURE 

QNETPP 

(W/SQ.M) 

0. 13282E+05 

RADIOSITY 

TEXTBOOK  [1] 
RADIOSITY 

1 

0.60000 

(KELVIN) 

1273.30 

(W/SQ.M) 

0.14016E+06 

(W/SQ.M) 

0.14003E+06 

2 

0.60000 

1273.30 

0.84339E+04 

0. 14339E+06 

0. 14326E+06 

3 

0.60000 

1273.30 

0. 15274E+05 

0. 13883E+06 

0. 13872E+06 

4 

0.60000 

1273.30 

0. 34994E+05 

0. 12569E+06 

0. 12557E+05 

5 

1.00000 

293.30 

- . 11659E+06 

0.41952E+03 

0.41952E+06 

SAMPLE  PROBLEM  (2)  -  CYLINDRICAL  CAVITY  WITH  SURFACE  1 
AT  1273.3  K,  WITH  SURFACES  2,3.  AND  4  INSULATED  (Q”=O.0)  , 
REJECTING  HEAT  TO  BLACKBODY  SINK  OF  293.3  K 


RADIATION  ANALYSIS 


SURFACE  EMI SSI VI TY  TEMPERATURE  QNETPP 

(KELVIN)  (W/SQ.M) 


TEXTBOOK  [1] 
RAD I OS I TY  RADI OS I TY 

(W/SQ.M)  (W/SQ.M) 


1  0.60000 
2  0.60000 

3  0.60000 

4  0.60000 

5  1.00000 


1273.30 
1093.27 
1004.92 
684 . 84 
293.30 


0.50442E+05 
O.OOOOOE+OO 
0 . OOOOOE+OO 
O.OOOOOE+OO 
- . 38822E+05 


0. 11539E+06 
0.80988E+05 
0.57815E+05 
O. 34751E+05 
0.41952E+03 


0. 11532E+05 
0.81019E+05 
0.57885E+05 
O. 34767E+05 
0.41952E+05 


k 


4 
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SAMPLE  PROBLEM  (3) 

(PROBLEM  16,  PACE  278  OF  REFERENCE  [3]) 

TWO  CONCENTRIC  CYLINDERS,  INNER  CYLINDER  AT  1800  K, 
ENVIRONMENT  AT  300  K 


THIS  PROBLEM  TESTS  THE  TWO-SIDED  SURFACE  OPTION  IN  SUBROUTINE  RADHT 
TEXT  BOOK  [3]  TEMPERATURE  OF  SURFACES  2,3  (OUTER  CYLINDER)  =  1048  K 


NUMBER  OF  ITERATIONS  =  15 

(RELATIVE  ERROR  =  1.0000E-03) 

RADIATION  ANALYSIS  FOR  ALL  SURFACES 


SURFACE  EMISSIVITY  TEMPERATURE  QNETPP  RADI  OS  I TY  IRRADIATION 

(KELVIN)  (W/SQ.M)  (W/SQ.M)  (W/SQ.M) 


1  0.22000  1800 

2  0.50000  1047 

3  0.17000  1047 

4  1.00000  300 


00 

1 . 1506E+05 

1 

49 

-1 . 1565E+04 

7 

92 

1 . 1544E+04 

1 

00 

0 .  OOOOE+OO 

4 

8694E+05  7 . 1880E+04 
9816E+04  9 . 1382E+04 
2003E+04  4.5919E+02 
5919E+02  4.5919E+02 


/  / 


<1  =  0.22  — '*  / — «z«0.50 

T|  =  1800  K 
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APPENDIX  B  -  COMPUTER  PROGRAM  LISTING 


PROGRAM  TO  DETERMINE  THE  SIZE  OF  A  MAIN  RADIATOR  AND 
A  SECONDARY  RADIATOR  ALONG  WITH  THE  MAXIMUM  OPERATING 
TEMPERATURES  OF  THE  PCU  AND  PAYLOAD 

WRITTEN  BY  DEAN  DOBRANICH,  NOV.  1986 

OUTPUT :  U[Jn1t^6  (MAIN).  UNIT  7  (VIEW  FACTORS).  UNIT  8  (ITERATIONS) 

COMMON/X/XCPCU , XCPL 
COMMON/EM  I  S/EM  I S ( 1 60 ) 

/lit  /ounDD  /  i  on  ) 


COMMON/'HT/ONETPPlieSi  ,T(100)  .  R£  1 00)  .  G(  1 00 ) 
COMMON/OUT /PMR .PSR.PPCUS.f  A"' ' 

DC  D 1  DCCO  APFfil 


130  ,Rl  iOW  _ 

HMK  , ran . rruua ,  PPCU  ,  PPLS , PPL ,  PCUWH  ,  PLWH  ,  REFPOWE  . 

PSR 1 , PSR2 . ARE AMR , AREASR , TAVGPCUS , TAVGPLS , COP . 

2  T I SMR . TOSMR . T I SSR , TOSSR , HAPCU , HAPL , CPCU , CPL , 

3  THPOWER  ,  THPNEW  _  ,  „„ 

COMMON/D  I MS/DPCUS , HPCUS . DPCU . HPCU . DPLS , HPLS , DPL , HPL . 

1  SD1 .SD2.SD3 

COMMON/AREA/APCUSL . APCUSC , APCUSR . APCUL , APCUC , APCUR . 

1  APLSL.APLSC.APLSR.APLL.APLC.APLR. 

2  APL . APCU . APLS , APCUS 

DIMENSION  IDSURF(100). IDCON(100), IST(20) 

DIMENSION  TSMR (50),TSSR(50) .TNODEMR(50) . T NODE SR (50 ) 

DIMENSION  RADMR(5M , R AD SR (50) . POAMR (50 ) ,POASR(50) .AMR (50) ,ASR(50) 
DIMENSION  RNOOEMR ( 50 ) ,RNOD£SR(50) 

DIMENSION  TPCU (4) . TPL (4 ) , RPCU( 4 ) . RPL (4 ) , FPCU( « , A) . FPL( A . 4 ) 
DIMENSION  ECPCU(A) ,ECPL(4) 


-3.141592654 


READ(5 , • )NMR 
RE AD ( 5 , • ) NSR 

READ(5 . • )EMR1 
REaD(5 . * 1EMR2 
READ(5.4)eSR1 

READ  5.»)ESR2  _  _ 

READ  5,*)EPCUS1L.EPCUS1C,EPCUS1R 
READ  5 . . )EPCUS2L .EPCUS2C .EPCUS2R 
READ  5..  EPLS1L.EPLS1C.EPLS1R 
REAO  5..)EPLS2L,EPLS2C,EPLS2R 
READ(5.4)EPCU 
READ(5.*)EPL 


C 


C 


C 


c 


c 


c 


c 


c 


c 

c 


READ(5.*)TSPACE 

READ(5. »)T I SMR 
READ(5 . * ) TOSMR 
READ(5.*)TISSR 
RE  AD  ( 5  .  •  ) TOSSR 

READ(5.4)RADMRI 

READ(5 . 4 1RADSR I 
RlAD(5.4)SD1 
READ ( 5 . • ) SD2 
READ ( 5 . • ) SD3 

READ(5 . • )DPCUS 
RE AD ( 5 . * IHPCUS 
READ(5.4)DPLS 
READ ( 5 . • ) HPLS 
RE AD ( 5 . • IDPCU 
READ(5 . • IHPCU 
READ(5 . * IDPL 
READ(5 . • )HPL 


READ(5. 
RE  AD (5 . 
RE  AD  ( 5  . 


JTHPOWER 

IEFFIC 

IFCCOP 


READ(5.4)FEPPCU 
READ(5.4)FPCUWHAC 
READC5 . « IFEPPL 
REA0(5 . * )FPLWHAC 

READ(5.4)DELTSR 
READ  5.4 )CPCU 
READ ( 5 , •  )CPL 


I DEBUG-0 

I F(NOUTERMAX . LT . 0) I  DEBUG- 1 
noutermax-abs(noutermax) 

NCV-2 

TDROPSR-ABS(T I SSR-TOSSR) 

I F [raDMR^ ' LE 1 DPCUS40 . 51RADMR 1-1 . 01 4DPCUS*0 . 5 

IF  ^ktDG^:T^-g^ELE1o5^S?^LO?D:D^S;?T?VE  COOLING  OF  PCU  ALLOWED 
I  F  (  FCCOP  .  NE  .  0 . 0  )  FPCUWHAC-0 . 0 
COP-0 . 0 
THPNEW-0 . 0 
REFPOWE-00 


COMPONENT  POWERS  (BEFORE  ADJUSTMENT  FOR  CONDUCTION  AND  CONVECTION) 
ELECPOW-THPOWER" E  F  F I C 
PMR-THPOWER-ELECPOW 
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PCUWH«FEPPCU*ELECPOW 
PSR1-FPCUWHAOPCUWH 
PP CUS-0.0 

ppcu-pcuwh-psri-ppcus 

PLWH-FEPPL  * ( E  LECPOW-PCUWH ) 

PSR2-FPLWHAC*PLWH 

PPLS-0.0 

PPL-PLWH-PSR2-PPLS 
PSR-PSR 1 +PSR2 
C 

C  COMPONENT  AREAS 

APCUSC-PI *DPCUS*HPCUS 
APLSC-PI *DPLS*HPLS 
APCUSL-PI *DPCUS*DPCUS*0 .25 
APLSL-PI *DPLS*DPLS*0.25 
APCUSR-APCUSL 
APLSR=APLSL 
APCUC-PI *DPCU*HPCU 
APLC-PI *DPL^HPL 
APCUL-P  I  *OPCU*DPCU*0 . 25 
APLL®P I *DPL*DPL*0 -  25 
APCUR*=APCUL 
APLR-APLL 
C 

APCUS-APCUSC+APCUSL+APCUSR 

APLS-APLSC+APLSL+APLSR 

APCU-APCUC+APCUL+APCUR 

APL-APLC+APLL+APLR 

C 

C  SET  UP  POINTERS 
I  ST ( 1 )-1 

1ST  i2)“=IST  ( 1  )+NMR 
I  ST ( 3 )- 1  ST ( 2  J+NMR 
IST(4)-IST(3)+NSR 
IST(5)-IST(4)+NSR 
IST(6)-IST(5)+1 
IST(7)-IST(6)+1 
IST(8)-IST(7)+1 
l$T(9)=IST(8)+1 
I  ST ( 1 0)« I  St (9)  +  1 
I  ST ( 1 1 )*IST( 10)  +  1 
I STi 12 )*l ST ( 1 1 )+1 
!$T(13)«IST(12)+1 
JST(14)-IST(13)+1 
I  ST ( 1 5 )- 1  ST  M  4)  +  1 
IST(16)-IST(15)+1 
I  ST ( 1 7 )» I $T i 1 6  )  +  1 
I  ST ( 1 8 )— I  ST ( 1 7 )  +  1 
IST(19)«=IST(18)  +  1 
C  TOTAL  NUMBER  OF  SURFACES 
NS® I  ST ( 1  9  ) 


C  SET  EMIS,  IDSURF.  AND  IDCON  ARRAYS 
C 

DO  5  K=1 , NS 
R  ( K )  - 1  . 0 
5  IDCON(K)-0 
C 

C  MAIN  RADIATOR 
KK-IST(2) 

DO  10  K— I  ST ( 1 ) .  1ST ( 2 )  —  1 
EMIS(K)-EMR1 
EM  I S (KK )-EMR2 
lDSURF(K)-0 
I DSURF (KK )— 0 
I  DCON(K)«KK 
I DCON ( KK )— K 
KK-KK+1 
10  CONT INUE 
C 

C  SECONDARY  RADIATOR 
KK- I  ST ( 4 ) 

DO  15  K— I  ST ( 3 ) . 1ST ( 4 )  —  1 
EMI S(K )— ESR1 
EM  I S (KK )*E  SR2 
I DSURF (K )-0 
I DSURF ( KK )=0 
I DCONfK KK 
I  DCON(KK)— K 
KK-KK+1 
15  CONT INUE 
C 

C  PCU  SHELL  LEFT  END 
KK— 1ST ( 6 ) 

DO  20  K— I  ST ( 5 ) .  I  ST ( 6 )  — 1 
EMIS(K)-EPCU$1L 
EMIS(KK)=EPCUS2L 
I DSURF ( K )— 2 
I DSURF ( KK )— 2 
I DCONfK )— KK 
I DCON ( KK — K 

ONETPP(K)— PPCUS/APCUSL/2 . 0 
ONETPP(KK)-ONETPP(K) 

KK-KK+1 
20  CONTINUE 
C 

C  PCU  SHELL  CURVED  SURFACE 
KK- 1  ST ( 8 ) 

DO  25  K— I  ST ( 7 ) , I  ST ( 8 )  — 1 
EM  I S (K )— EPCUS1 C 
EM  I S (KK )— EPCUS2C 
I DSURF(K)-2 
I DSURF (kk)*2 
I DCON ( K )— KK 


IDCON(KK)  —  K 

QNETPP(K)-PPCUS/APCUSC/2 .0 
QNETPP(KK)«QNETPP(K) 

KK-KK+1 
25  CONTINUE 
C 

C  PCU  SHELL  RIGHT  END 
KK=IST( 10) 

DO  30  K«IST(9) , t  ST ( 10)  —  1 

EMIS(K)«EPCUS1R 

EMIS(KK)«EPCUS2R 

IDSURF(k)-2 

IDSURF(KK)-2 

IDC0N(K)-KK 

IDCON(KK)— K 

GNETPP(K)«PPCUS/APCUSR/2 .0 
ONETPP (KK  ) —ONE TPP (K ) 

KK-KK+1 
30  CONTINUE 
C 

C  PAYLOAD  SHELL  LEFT  END 
KK«IST(12) 

DO  35  K«IST(1 1) . I  ST ( 1 2 )  —  1 

EMIS(K)-EPLS1L 

EM  I S(KK)-EPCUS2L 

I  DSURF (K)-2 

IDSURF(KK)-2 

IDC0N(K)-KK 

I DCON(KK)— K 

QNETPP(K)-PPLS/APLSL/2 . 0 
ONETPP (KK)-QNETPP(k') 

KK-KK+1 
35  CONTINUE 
C 

C  PAYLOAD  SHELL  CURVED  SURFACE 
KK=IST( 14) 

DO  40  K«IST(13) . IST(14)-1 

EM  I S (K )-EPLS1 C 

EMIS(KK)-EPCUS2C 

ID$URF(K)-2 

I DSURF (KK)*2 

IDCON(K)-KK 

I DCON(KK)=— K 

ONE  TPP ( K )«PPLS/APLSC/2 . 0 
ONE  TPP  ( KK  )»=ONETPP  ( K ) 

KK-KK+1 
40  CONTINUE 
C 

C  PAYLOAD  SHELL  RIGHT  END 
KK«=IST(  16) 

DO  45  K- I  ST ( 1 5) . I  ST ( 16 )-1 
EMI  S(K)«=EPLS1R 


C 

C 


C 

C 


C 

C 


C 

C 


c 


c 

c 

c 

c 


EMIS(KK)«EPCUS2R 
I DSURF ( K )=2 
I DSURF ( KK )=2 
I DCON ( K )— KK 
I  DCON(KK)— K 

ONE  TPP  ( K  )«=PPLS/ APLSR/2 . 0 
ONETPP (KK )=ONETPP (K) 
KK-KK+1 
45  CONTINUE 


PCU 

EM  I S ( ! ST (17) )-EPCU 

I DSURF ( I  ST ( 1 7 ) )-1 

ONE TPP ( IST(17))=PPCU/APCU 

PAYLOAD 

EM  I S ( I  ST ( 1 8 ) )=EPL 
I  DSURF 7  I  ST ( 1 o ) )  — 1 
ONETPP ( IST(18))“PPL/APL 

SPACE 

EM  I S ( I  ST ( 19) )«1  . 0 
I DSURF ( I  ST  C 1 9) )— 0 
T( I  ST ( 19) )  — T  SPACE 


INITIAL  GUESS  FOR  EFFECTIVE  TEMPERATURES 
DO  60  J-1 , NMR 
60  TSMR( J )— TSPACE 
DO  65  J-1 , NSR 
65  TSSR(J)—T SPACE 


RADMR( 1 )-RADMR I 
RADSR( 1 )=RADSRI 
PNMR-PMR/NMR 
PNSR—PSR/NSR 
EAVGMR- ( EMR 1 +EMR2 )/2 . 0 
EAVGSR-(ESR1+ESR2)/2 .0 
WRITE  SURFACE  ID  ARRAYS 
CALL  I DOUT ( NS . NMR , NSR , I 
DETERMINE  VIEW  FACTORS  FOR 
THE  SIZE  OF  THE  MAIN  AND 
CALL  VIEW1 (NS, 1ST) 

SET  UP  ARRAYS  FOR  SOLVING 
DO  120  1-1,4 
RPCU( I )-1 .0 
120  RPL ( I )-1 .0 

ECPCUt 1 )-EPCUS1 L 
ECPCU ( 2 )-EPCUS 1 C 
ECPCU ( 3 )— EPCUS 1 R 
ECPCU ( 4 )— EPCU 
ECPLfl )-EPLS1L 
ECPL(2)-EPLS1C 


DSURF, I DCON) 

SURFACES  NOT  DEPENDENT  ON 
SECONDARY  RADIATORS 

COMBINED  RADIATION/CONDUCT 


ON 
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ECPL(3)-EPLS1R 

ECPL(4)-EPL 

CALL  FASSIGN(0, I  ST , FPCU) 
CALL  FASSIGN(l  ,  1ST, FPL) 


DO  600  ICV-1.NCV41 
00500  IOUTER-1 .NOUTERMAX 

FIRST  ESTIMATE  THE  SIZE  OF  THE  MAIN  RADIATOR 

CALL  RADS  I ZE ( NMR , PMR , T I SMR , TOSMR , EAVGMR , T SMR , RADMR .AMR , 

1  TNODEMR .RNODEMR) 

ASSIGN  TNODEMR  ARRAY  TO  T  ARRAY 
KK- I  ST (2) 

Do' 75  J-IST(I), I ST(2 )— 1 
T(J)«=TNODEMR(K) 

T  ( KK ) =TNODEMR ( K ) 

K-K+1 
KK-KK+1 
75  CONTINUE 

NOW  ESTIMATE  SIZE  OF  THE  SECONDARY  RADIATOR 

CALL  RADS  I ZE (NSR , PSR . T ISSR , TOSSR . EAVGSR . TSSR , RADSR , ASR , 

1  TNODESR .RNODESR) 

ASSIGN  TNODESR  ARRAY  TO  T  ARRAY 
KK* I  ST ( 4 ) 

Do's0  J-ISTE3) . IST(4)-1 
T(J)*TNODESR(K) 

T ( KK )=TNODESR ( K ) 

K=K41 
KK*KK4l 
80  CONTINUE 

NOW  DETERMINE  VIEW  FACTORS  FOR  ALL  SURFACES  BASED  ON  EST I  MATED 
AREAS.  THEN.  SOLVE  FOR  RADIOSITIES,  ONETPP  S.  AND  T  S. 

CALL  VIEW2(NS , NMR .NSR .  I  ST , RADMR . AMR . RADSR , ASR ) 

NITER-NRHMAX 

CALL  RADHT (NS , IDSURF. I DCON , CONVRH .NITER) 

DETERMINE  POA  BASED  ON  CALCULATED  ONETPP'S  OF  MA I N  AND  SECONDARY  RADIATORS 
WHERE:  POA  -  ONETPP(SIDE  1)  4  ONETPP(SIOE  2) 

ALSO.  CHECK  FOR  CONVERGENCE  OF  ONETPP ‘S 

MAIN  RADIATOR  FIRST 
RELERRM-0.0 
POWT OTMR*=0  0 
KK* I  ST ( 2) 

K=1 


C 

C 


C 


C 


C 

C 


c 

C 


c 

C 

C 

c 


DO  100  J-IST(t) . I  ST (2)-l 
POAMR (K)-0NETPP(J)40NETPP(KK) 
PNMRNEW-POAMR ( K ) • AMR ( K ) 
POWTOTMR-POWTOTMR+PNMRNEW 
RELERR-ABS(PNMR-PNMRNEW)/PNMR 
RE  LERRM=AMAX 1 (RELERR .RELERRM) 
K*K41 
KK=KK41 
00  CONTINUE 


SECONDARY  RADIATOR  LAST 
POWTOTSR-0 . 0 
KK*  I  ST  (  4  ) 

K=  1 

DO  105  J= I  ST (3) , I  ST ( 4 )-1 

P0ASP(K)-0NETPP(J)40NETPP(KK) 
PNSRNEW-POASR(K) .ASR(K) 

POWT  OT  SR*POWT  0TSR4PNSRNEW 
RELERR=ABS (PNSR-PNSRNEW)/PNSR 
RELERRM=AMAX1 (RELERR .RELERRM) 
K-K41 
KK-KK41 
105  CONT INUE 


WRITE(8. 88) IOUTER. NITER. RELERRM, POWTOTMR , POWTOTSR 
8e  FORMATC  ’ ,/,4X, I2.3X, I2.3X, 1PE12.4 .3X.E12.4.3X.E12.4) 

I  F  (RELERRM  .  LT  .  CONVOUT  )G0  TO  505 


1 


■( I  DEBUG . EO- 1  .AND .  I CV . GT . NCV ) THEN 

4LL  ROUT  JnM^  .  NSR  .  TNODEMR  ,  TNODESR  ,  TSMR  .  TSSR  , 

POAMR  .  POASR  .  RADMR  ,  RADSR  .  AMR  .  ASR  .  RNODEMR  ,  RNODESR  ) 


CALL  OUT  1 (NS) 

ELSE 

END  IF 


CALCULATE  NEW  EFFECTIVE  TEMPERATURES  (TSNEW’S) 
CALL  TSNEW( NMR .EAVGMR . TNODEMR . POAMR , TSMR) 
CALL  TSNEW( NSR .EAVGSR , TNODESR . POASR .TSSR) 
AND  REPEAT  THE  ENTIRE  PROCESS... 


IF  CONVERGENCE  NOT  MET 


500  CONTINUE 

WR I T  E ( 6 , 499) 

499  FORMAT ( '  •,/.4X.,-*« 
1  NOT  CONVERGE 


.  WARNING  -  OUTER 

/) 


ITERATION  DID 


505  CONTINUE 


ARE AMR-P I  * (RADMR ( NMR4 1 
AREASR-P I  * ( RADSR (NSR41 


*  *2-RADMR 1**2) 

*  *2-RADSR I »*2) 


r 
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NOW  ACCOUNT  FOR  CONDUCTION  BETWEEN  PCU  AND  PCU  SHELL 

USE  AVERAGE  TEMPERATURE  OF  SHELL  TO  ALLOW  FOR  MAXIMUM  POSSIBLE 
TEMPERATURE  REDUCTION  DUE  TO  CONDUCTION 
TAVGPCUS-(T( IST(5))*APCUSL+T( I  ST ( 7 ) ) *APCUSC+ 

1  T  ( I  ST ( 9 ) ) * APCUSR )/APCUS 

TAVGPLS-(T( I  ST ( 1 1  ) )*APLSL+T  < I  ST (13)) .APLSC+ 

1  t(IST(15))*APLSR)/APLS 

TPCU(1)-TAVGPCUS 
TPCU(2)=TAVGPCUS 
TPCU(3)-TAVGPCUS 
TPL  ( 1  )  -T  AVGPLS 
TPL(2)-TAVGPLS 
TPL(3)— TAVGPLS 

WRITE  FINAL  OUTPUT  FOR  NO-SHELL  CONDUCTION  CASE 
CALL  OUT2( 1ST) 

USING  MODIFIED  REGULA  FALSI  METHOD  TO  FIND  ROOT  OF  FUNCTION  FUN 
SUCH  THAT  FUN (ROOT ) -WASTE  HEAT-RADIATED  HEAT-CONDUCTED  HEAT-0 
FIRST  PCU  (HI  AND  H2  PROVIDE  AN  INTERVAL  CONTAINING  A  ROOT) 
IF(CPCU.LE.0.0)GO  TO  150 
Hl-0 . 0 
H2-PPCU 
ROOTN-PPCU 

CALL  FUN ( PPCU  .  APCU  .  ECPCU  .  RPCU  .  FPCU  .  TPCU  .  CPCU , XC . H 1 , FH | 

CALL  FUN ( PPCU , APCU .ECPCU . RPCU . FPCU , TPCU , CPCU , XC , H2 . CH) 

ROOT  0—H2 

DO  125  I M—1 .20 

GMF-GH-FH 

I F(ABS(GMF) .LT. 1 ,0E-15)GO  TO  150 
ROOTN=(GH*H1— FH*H2)/GMF 
CHECK  FOR  ERROR 

I F ( ABS( (ROOTN-ROOTO)/ROOTN) . GT . 0 . 002)GO  TO  140  . 

CALL  FUN (PPCU . APCU , ECPCU , RPCU . FPCU , TPCU . CPCU . XC , ROOTN , FROOT ) 

1 40  CALL°FUN ( PPCU , APCU , ECPCU . RPCU . FPCU . TPCU , CPCU . XC , ROOTN . FROOT ) 
I F ( FH* FROOT . GT . 0 . 0) GO  TO  130 
H2-R00TN 
GH-FROOT 
FH-0.5*FH 
GO  TO  135 
130  HI -ROOTN 
FH-FROOT 
GH-0. 5*GH 
135  CONTINUE 

ROOTO-ROOTN 
125  CONTINUE 
GO  TO  170 


150  CONTINUE 

I F ( CPCU . EO . 0 . 0) GO  TO  175 
TPCU  (  4  )=AMAX  1  ( TPCU (  1  )  ,TPCU(2)  . TPCU ( 3 ) )-CPCU 

CALL  FUN (PPCU .APCU . ECPCU . RPCU . FPCU . TPCU , CPCU . XCPCU . ROOTN . FROOT ) 
170  PPCU=ROOTN 

PPCUS=PCUWH— PSR 1 —PPCU 
ONETPP(IST(17)) —PPCU/ APCU 
ONE T PP (  I  ST (5) )=PPCUS/APCUS 
ONETPP( IST(7) )-PPCUS/APCUS 
ONE TPP ( I  ST ( 9  j )=PPCUS/ APCUS 
ONETPP ( IST(6))«0.0 
ONETPP ( I  ST ( 8 ) )-0 • 0 
ONETPP ( I  ST ( 10) )-0 .0 
175  CONTINUE 

REPEAT  FOR  PAYLOAD 

I F ( CPL . LE . 0 . 0 ) GO  TO  250 
HI  =0.0 
H2-PPL 
ROOTN-PPL 

CALL  FUN(PPL.APL.ECPL,RPL.FPL.TPL.CPL.XC.H1 ,FH) 

CALL  FUN(PPL.APL.ECPL,RPL.FPL.TPL,CPL,XC.H2,GH) 

ROOT  0=H2 
DO  225  I M—1 ,20 
GMF-GH-FH 

I F ( ABS (GMF )  .  LT  .  1  . 0E-1 5 ) GO  TO  250 
R00TN-(GH*H1-FH*H2)/GMF 

CHECK  FOR  ERROR  v  ^ 

I F ( ABS ( (ROOT N-ROOTO) /ROOTN) . GT . 0 . 002) GO  TO  240 
CALL  FUN ( PPL, APL.ECPL.RPL. FPL. TPL, CPL. XC. ROOTN, FROOT) 

240  CALL°FUN?PPL,APL ,ECPL ,RPL , FPL, TPL , CPL. XC, ROOTN . FROOT) 

I  F ( FH* FROOT . GT . 0 . 0 ) GO  TO  230 
H2-ROOTN 
GH-FROOT 
FH— 0 . 5*  FH 
GO  TO  235 
230  H 1 —ROOTN 
FH-FROOT 
GH— 0 . 5*  GH 
235  CONTINUE 

ROOTO-ROOTN 
225  CONTINUE 
GO  TO  270 

250  CONTINUE 

I F ( CPL . EQ . 0 . 0 ) GO  TO  275 

TPL ( 4 )— AMAX 1(TPL(1) ,TPL(2),TPL(3) )— CPL  , 

CALL  PUN( PPL, APL.ECPL.RPL. FPL, TPL. CPL, XCPL, ROOTN, FROOT) 

270  PPCU— ROOTN 
PPL-ROOTN 
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c 

c 

c 

c 


c 


c 

c 

c 


c 

c 


c 


PPLS— PLWH— PSR2— PPL 
ONETPP ( J  ST ( 1 8 ) )-PPL/APL 
QNETPPf IST(11))«PPLS/APLS 
ONETPPi  i ST  I  1 3 ) J-PPLS/APLS 
ONETPPI I  ST ( 15))-PPLS/APLS 
ONETPP ( IST(12))«0.0 
ONETPP l IST(14))-0.0 
ONETPP ( IST(16))-0.0 
275  CONTINUE' 

RE-SOLVE  RADIATION  HEAT  TRANSFER  FOR  ENTIRE  PLATFORM  USING 

PCU  AND  PL  WASTE  HEAT  VALUES  CORRECTED  FOR  CONDUCTION  TO  SHELL 

NITER-NRHMAX 

CALL  RADHT (NS , IDSURF, I DCON , CONVRH , N I  TER) 

TAVGPCUS-(T( I  ST (5) ) *APCUSL+T ( I  ST ( 7 ) ) * APCUSC+ 

1  T  ( I  ST ( 9 ) ) * APCUSR)/APCUS 

TAVGPLS-(T(  I  ST  ( 1  1  ) )*APLSL+T(  I  ST ( 1 3) ) * APLSC+ 

1  T  ( I $T ( 15) ) *APLSR)/APLS 

ADJUST  MAX  SEC  RADIATOR  TEMP  SO  THAT  IT  DOES  NOT  EXCEED  MAX  OF 

PL  OR  PCU  TEMPS  (OTHERWISE,  A  REFRIGERATION  UNIT  WOULD  BE  NEEDED) 
I F (NCV . EO . 0 . OR . r CCOP . NE . 0 . 0)GO  TO  599 
TM1-T( I  ST ( 1 7 ) ) 

I F (PSR1 .  LT  .  1 . 0) TM1— 0 . 0 
TM2«T( ! ST ( 18) ) 

I F ( PSR2 . LT . 1 . 0) TM2-0 . 0 

I F( T I SSR  GT . TOSSR ) THEN 

T I SSR=AMAX 1  (  TM1 , TM2)-DELTSR*0 . 97 

TOSSR-T  I  SSR-TDROPSR 

ELSE 

TOSSR-AMAX1 (TM1 , TM2)-DELTSR *0 . 97 
T ISSR-TOSSR-TDROPSR 
END  I  F 

CALCULATE  REOU I  RED  HA  FOR  CONVECTION  COOLING  OF  PCU  AND  PL 
(IF  HA  IS  NEGATIVE,  A  REFRIGERATION  UNIT  IS  REQUIRED) 

599  TCAVG=(T I SSR+ TOSSR )/2 . 0 

HAPCU=PSR 1 /( T ( I ST( 1 7 ) )-T  CAVG) 

HAPL«PSR2/(T( l$T(18) j-TCAVG) 

IF  A  REFRIGERATOR  IS  USED  TO  COOL  PAYLOAD: 

IFfFCCOP.EO. 0.0. OR. I CV . GT . NCV) GO  TO  600 

COP-FCCOP* ( 1 • 0/( ( T I SSR/( T ( I  ST ( 1 8) )-0 . 96*DELTSR ) )-l . 0) ) 

REFPOWE=PSR2/COP 

REFPOWT-REFPOWE/EFFIC 

THPNEW=THPOWER+REFPOWT 

ELECPOW-THPNEW*EFFIC 

PMR=THPNEW— ELECPOW 

PNMR=PMR/NMR 

PSR 1 =0 . 0 

PPCUS=O.0 

PCUWH=FEPPCU* ELECPOW 


C 

C 


C 


C 

C 

C 

C 

C 

C 

C 


PPCU-PCUWH 

ONETPP ( 1ST (17))*PPCU/APCU 
ONETPP( I  ST ( 5 ) )=PPCUS/APCUS 
ONE T PP (  I  ST  (7)  )«=PPCU$/APCUS 
ONETPPi  1ST  ( 9 )  )*=PPCUS/APCUS 
ONETPPi IST(6))=0.0 
ONETPP(  1ST  C  8 ) )=0 . 0 
ONE TPP ( I  ST ( 1 0 ) )=0 . 0 
600  CONTINUE 

NOW  WRITE  ALL  FINAL  OUTPUT 
CALL  FOUT(NS) 

CALL  ROUT (NMR ,NSR , TNODEMR , TNODESR , TSMR , TSSR . 

1  POAMR , POASR , RADMR , RADSR , AMR , ASR , RNODEMR , RNODESR ) 

CALL  OUT  1 (NS) 

CALL  OUT2 ( 1ST) 

STOP 

END 


* - SUBROUTINES - * 

SUBROUT  I NE  RADS  I ZE ( NN , P , T I S , TOS , EAVG ,TS. RAD . A , TNODE , RNODE ) 

THIS  SUBROUTINE  DETERMINES  THE  REOU I  RED  SIZE  OF  THE  RADIATOR  NODES 
BASED  ON  AN  ESTIMATED  EFFECTIVE  BACKGROUND  TEMPERATURE  ( TS ) - 

DIMENSION  RAD ( NN+1 ) ,TS(NN) , TNODE (NN) . RNODE (NN) , A ( NN ) 

PI-3. 141592654 


C 

C 

C 

C 


c 

c 

c 

c 

c 


c 


I F (ABS(T I S-TOS) . GT .2.0) THEN 
FOR  A  BRAYTON  CYCLE: 


C2— 1  .  0/4 . 0/P  I  /S  I  GMA/E  AVG 
C2—  1  .40373E06/EAVG 
C1—C2  *P/ABS  f T I S-TOS ) 

LET  EACH  NODE  HAVE  THE  SAME  TEMPERATURE  DROP;  DETERMINE 
NODE  RADII  SUCH  THAT  POWER  DEPOSITED  IN  NODE  EQUALS  POWER 
REMOVED  FROM  BOTH  SIDES  OF  NODE  VIA  RADIATION 
DT=(TOS-T IS)/NN 
FOR  OUTS  IDE- IN  FLOW 
S I GN=1 .0 


FOR  INSIDE-OUT  FLOW 

IF(DT.LT.0.0)SIGN— 1 .0 
DO  5  1-1 ,NN 
T0-TIS+DT*( 1-1) 

T 1 =T I S+DT* I 
ARG-AMIN1 (T0.T1 ) 

TS  CAN  NOT  BE  GREATER  THAN  ARG 


I F ( TS ( I ) . GE . ARG ) TS ( I 
TRAT IO=(TS(!)+T1)/(T_ ,  ,  ,  ,  , 

RAD( 1+1 )  — ( RAD ( I ) *  *  2+S I GN*C1 *0.5/tS( I ) *  *3* 


I )— ARG*  0 . 998 
TS( I ) — T 1 )*(TS(I )-T0)/(TS( I )+T0) 


:  44 


c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 


c 

c 

c 

c 

c 

c 

c 

c 


( ALOG( TRAT I 0)+2 . 0* AT  AN ( T1/TS( I ))- 
2.0*ATAN(T0/TS( I )))v  *  " 


5  CONTINUE 

DETERMINE  AVERAGE  NODE  TEMPERATURE  BASED  ON  TEMPERATURES  AT  END 
OF  EACH  NODE  AND  THE  NODE  EFFECTIVE  BACKGROUND  TEMPERATURE; 
ALSO,  CALCULATE  R  AT  WHICH  AVERAGE  TEMPERATURE  OCCURS 

C3-C1 *P I 
ADT-ABS(DT) 

DO  10  1-1 ,NN 

A (  l^-Pi *(RAD( 1+1 )**2-RAD( I )**2) 


DE ( I )-(TS( I )**4— 2 . 0*C3*ADT/A( I  )  )  *  *0 . 25 
T I S+DT  * ( 1-1 ) 


TNOi 
T0. 

T 1-TNODE ( I 
TRAT IO=(T 
RNODE ( I )- 

1 

2 

10  CONTINUE 


S( 1 )+T1)/(T$( 

(RAD( 

) *  *2+S 1 

(ALOG 

TRAT 10) 

2.0*An 

Un(T0/T 

S(  I  ) )~ 


RETURN 

ELSE 

FOR  A  RANKINE  CYCLE: 

PONN-P/NN 

T  I  S4-T I S*  *4 

C4-2 . 0*5 . 669E~8*EAVG 

DO  15  1-1 , NN 

A ( I )«P0NN/C4/(TIS4-T$( I )**4) 

RAD( 1  +  1 )«(A( I )/PI+RAD( I ) *  *2 ) *  *0 . 
RNODE ( I )-0. 5* ( RAD ( I )+RA&( 1+1)) 
TNODE ( I j-TlS 
15  CONTINUE 


END  I  F 
RETURN 


SUBROUTINE  T SNEW( NN , EAVG , TNODE , POA , TS) 

THIS  SUBROUTINE  FINDS  NEW  EFFECTIVE  BACKGROUND  TEMPERATURES 
DIMENSION  TNODE (NN) ,POA(NN) ,TS(NN) 

C1-1/2/S I GMA/EAVG 
C1-8 . 820E06/EAVG 

WEIGHTING  FACTOR  USED  TO  DAMPEN  OSSCILATIONS  AND  SPEED  CONVERGENCE 
WF-0 . 8 


C 


c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 


DO  5  1-1 , NN 

TS( I )-WF ♦ ( TNODE ( I ) • *  4-POA ( I ) *C1 ) * *0 . 25+ ( 1 . 0-WF)*TS( I ) 
5  CONTINUE 

RETURN 

END 


SUBROUTINE  RADHT ( N ,  IDSURF,  I DCON . CONV . N I  TER ) 
COMMON/CV I EW/F (100,100) 

COMMON/EM  I  S/EM  I S ( 1 00 ) 

COMMON/HT/GNETPP ( 100) , T ( 1 00) . R( 100) , G( 1 00 ) 
DIMENSION  V( 1 80 ) , IDSURF(N), IDCON(n) 

CALCULATE  RAD  I  OS  I T I ES ,  I RRAD I  AT  I ONS , AND  HEAT  FLUXES 
USING  NETWORK  METHOD  FOR  RADIATION  HEAT  TRANSFER 
(REF.  HOLMAN,  HEAT  TRANSFER) 


IDSURF  -  0  FOR  SURFACES  WITH  KNOWN  TEMPERATURE 

IDSURF  -  1  FOR  SURFACES  WITH  KNOWN  HEAT  FLUX 

IDSURF  =  2  FOR  TWO-SIDED  SURFACES  WITH  T1-T2,  QT-01+02 


I DCON  ARRAY  DEFINES  WHICH  TWO-SIDED  SURFACES  ARE  CONNCECTED ; 
SIDE  2  OF  A  CONNECTION  MUST  BE  INDICATED  AS  A  NEGATIVE. 
EXAMPLE:  IF  SURFACE  5  AND  22  ARE  CONNCECTED  - 

I DCON ( 5)-22  SIDE  1 

IDCON(22)— 5  SIDE  2 

SIGMA-5.669E-8 

USING  GAUSS-SE I  DEL  TO  SOLVE  NXN  EQUATIONS 
INITIAL  GUESS  FOR  RADIOS  I  TIES 
DO  51  1-1 , N 
51  V( I )— R ( I ) 

DO  99  IT-1 .NITER 

CALCULATE  R  FOR  SURFACES  WITH  KNOWN  TEMPERATURE 
DO  61  1-1 ,N 

I F ( I DSURF ( I ) . NE . 0 ) GO  TO  61 
X— 1 . 0— EM  I S ( I ) 

XX-1 . 0-F ( I , l )*X 

SUM-0.0 

DO  71  J— 1 . N 

!F(J .EO. I )GO  TO  71 

SUM—SUM+F ( I , J)*R(J) 

71  CONTINUE 

R ( I )-( X+SUM+EM I S( l)*SIGMA*T( I )**4.0)/XX 
61  CONTINUE 

CALCULATE  R  FOR  SURFACES  WITH  KNOWN  HEAT  FLUX 
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ooooo 


DO  62  1-1 ,N 

IF( IDSURFf I ) . NE . 1 ) GO  TO  62 
Y— 1 .0/(1 . 0-F ( 1,1)) 

SUM-0 . 0 
DO  65  J-1.N 
I F(J .EQ. I )G0  TO  65 
SUM— SUM+F ( I , J ) *R ( J ) 

65  CONTINUE  ,  ,, 

R ( I )-Y.(SUM+ONETPP( I )) 

62  CONTINUE 

CALCULATE  R  FOR  SURFACES  THAT  ARE  OPPOSITE  SIDES  OF  THE  SAME 
STRUCTURE  SUCH  THAT  0( TOTAL)  -  0(1)  +  Q(K)  AND  T ( I )  -  T(K). 

K  -  IDCON(I)  DEFINES  SURFACE  CONNECTIONS 

DO  63  1-1 ,N 

IF( IDSURFf I ) .NE .2)G0  TO  63 

K- I DCON ( I  ) 

I  F(K . LE . 0)GO  TO  63 

o*(xi-i.0)) 

B 1 — 1 .0/(1 . 0-F ( K , K ) ) 

X2-1 ■ 0/EMI S(K) 

A2-1  0/(X2-F(K.K)»(X2-1 .0)) 

62-1 .0/(1  - 0-F ( 1,1)) 

SUM1-0 .0 

SUM2-0 . 0 

DO  10  J-1  ,N 

1 F ( J . EO . I )G0  TO  10  , 

SUM1— SUM1+F ( I , J ) *R ( J ) 

SUM2-SUM2+F(K,j)*R(j) 

10  CONTINUE 
SUM3-0.0 
DO  12  J-1 ,N 

I F ( J . EO . K )G0  TO  12  ,  ,  > 

SUM3-SUM3+ ( F ( I ,J)+F(K, J))*R(J) 

12  R?niA1  •(  (X1-1 . 0  )  *SUM1  —  (  X2— 1  .  0) »SUM2+B  1*X2*(0NETPP(  I  )  + 

1  0NETPP(K)+SUM3) )/(1 . 0+A1 *B1 *X2) 

SUM 1-0 . 0 

SUM2-0.0 

DO  15  J— 1 , N 

I F ( J . EO . K ) GO  TO  15 

SUM1 — SUM1+F (K,J)*R(J) 

SUM2— SUM2+F f I , j).R(j) 

15  CONTINUE 
SUM3-0.0 
DO  17  J-1 ,N 

I F ( J .EO. I )G0  TO  17  _  ,  , 

SUM3-SUM3+(F( I ,J)+F(K,J)1*R(J) 


17  R?K)-A2<((X2-1 . 0 ) * SUM1  —  ( X 1  —  1  . 0) . SUM2+B2 • X 1  * (ONETPP ( I )  + 
1  QNETPP(K)+SUM3))/(1 -0+A2«B2*X1) 

63  CONTINUE 


CHECK  FOR  CONVERGENCE 

DIFF-0.0 
DO  81  1-1 , N 
DN=0  0 

I F (ABS ( V( I )) .GT.0.0)DN*ABS((R( I )-V( I ))/V( I )) 

81  D I FF=AMAX 1 ( D I FF , DN ) 

I F ( D I FF . LT . CONV. AND . D I FF . NE . 0 . 0)GO  TO  101 
DO  82  1-1 ,N 

82  V ( I )-R ( I ) 

99  CONTINUE 

101  CONTINUE 


DETERMINE  IRRADIATIONS  G,  FOR  EACH  SURFACE 
DO  120  1-1 , N 
SUM-0.0 
DO  110  K— 1 .N 
110  SUM=SUM+F( I , K)*R(K) 

120  G( 1 )— SUM 

DETERMINE  ONETPP  AND  t  FOR  EACH  SURFACE  (WATTS/M^2,  AND  K) 

DO  125  1-1 ,N  .  . 

I F ( I DSURF ( I ) . NE . 1 }QNETPP( I )-R( I )-G( I ) 

I  F (  IDSURFf I ) . NE . 0) THEN 
X  X— 0  0 

I F ( EM  I S ( I ).NE.0.0)XX-(1 ,0-EMIS( I ) )/EM I S ( I ) .ONETPP ( I )/SIGMA 
T ( I )  — ( XX+R ( I )/S I GMA ) *  *0 . 25 
ELSE 
END  I  F 

125  CONTINUE 

N I  TER- I T 

RETURN 

END 


SUBROUT  I NE  FUN ( WH ,A,EMIS,R,F,T.C,XC,0,FV) 

FUN  -  WASTE  HEAT  -  HEAT  RADIATED  -  HEAT  CONDUCTED 

USED  TO  COUPLE  RADIATION  AND  CONDUCTION  HEAT  TRANSFER 

DIMENSION  EM !S(4) ,R(4) . ONETPP (4) , F ( 4 , 4) ,T(4) 
DIMENSION  G ( 4 ) ,V(4) 
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CALCULATE  RADIOSITIES, IRRADIAT IONS. AND  HEAT  FLUXES 
USING  NETWORK  METHOD  FOR  RADIATION  HEAT  TRANSFER 

CONDUCTION  ENERGY  N v 

QC(TT)-C/3 . 0* ( (TT-T ( 1 ) )+(TT— T (2) )+(TT— T (3) ) ) 

S I GMA-5 . 669E— 8 
N-4 

CONV-1 . 0E-3 

I F(C . LT . 0 . 0)THEN 

I  FLAG-1 

NC-N 

ELSE 

I  FLAG-0 

NC-N— 1 

END  IF 

USING  GAUSS-SE I  DEL  TO  SOLVE  4X4  EOUATIONS 
INITIAL  GUESS  FOR  RADIOSITIES 
DO  51  1-1 , N 
GNETPPf I )— 0.0 
51  V( I )-R( I ) 

IF( l FLAG . EQ . 0)ONETPP(N)— O/A 
DO  99  IT-1 ,15 

CALCULATE  R  FOR  SURFACES  WITH  KNOWN  TEMPERATURE 
DO  61  1-1 ,NC 
X— 1 . 0-EM I $( I ) 

XX-1 . 0— F ( I , I )*X 
SUM-0 . 0 
DO  71  J-1.N 
I F( J . EQ . I )GO  TO  71 
SUM— SUM+F ( I .  J)*R(J) 

CONTINUE 

R ( I |-(X*SUM+EM I S( I ) *S I GMA*T ( I )**4.0)/XX 


CONTINUE 


CALCULATE  R  FOR  SURFACES  WITH  KNOWN  HEAT  FLUX 
DO  62  I  — N , N- I  FLAG 
Y-1 .0/(1 .0-F( 1,1)) 

SUM-0 . 0 
DO  65  J-1 ,N 
I  F ( J .EQ. I ) GO  TO  65 
SUM-SUM+F ( I , J) *R(  J) 

65  CONTINUE 

R ( I ) — Y  * 7  SUM+QNE  T  PP ( I ) ) 

62  CONTINUE 

CHECK  FOR  CONVERGENCE 

DIFF-0.0 


DO  81  1-1 , N 
DN-0  0 

! F (ABS ( V ( I ) )  GT .0-0)DN-ABS( (R( I )-V( I ) )/V( I ) ) 

81  DIFF-AMAXI(DIFF.ON) 
IF(DIFF.LT.CONV.An6.DIFF.NE.0.0)GO  TO  101 
DO  82  1-1 ,N 

82  V(l)-R(l) 

99  CONTINUE 

101  CONTINUE 


DETERMINE  IRRADIATIONS  G.  FOR  EACH  SURFACE 
DO  120  1-1 , N 
SUM-0.0 
DO  110  K— 1 , N 
110  SUM— SUM+F ( I ,K)*R(K) 

120  G( I )-SUM 


DETERMINE  ONETPP  AND  T  FOR  EACH  SURFACE  (WATTS/M«*2,  AND  K) 

DO  125  1-1 , NC 
125  ONETPP(I)— R(l)-G(l) 

IF(IFIAG.EO.0)THEN 
I  — N 

?f7em?S( I ) .NE .0 . 0 ) XX— ( 1 .0-EMI S( I ) )/EMIS( I ),QNETPP( I )/S I GMA 
T ( I )-( XX+R ( I )/S I GMA ) *  *0 . 25 
DETERMINE  VALUE  OF  FUNCTION 
FV— WH-QNETPP ( N ) * A— OC ( T (N ) ) 

ELSE 

0— ONETPP (N) *A 

TMAX—AMAX 1 ( T ( 1 ) , T ( 2) . T (3) ) 

XC— (WH-O)/ ( T ( 4 )— TMAX ) 

FV— 0 . 0 
END  IF 


RETURN 

END 

- *  *  *  * - 

SUBROUTINE  VIEWI(N.IST) 

CALCULATE  VIEW  FACTORS  FOR  ALL  SURFACES  NOT  DEPENDENT  ON 
SIZE  OF  MAIN  AND  SECONDARY  RADIATORS 
DIMENSION  I  ST (20 ) 

COMMON/CV I EW/F (100,100)  „  „ 

COMMON/AREA/APCUSL  ,  APCUSC  ,  APCUSR  ,  APCUL  .  APCUC  ,  APCUR  , 

1  APLSL . APLSC , APLSR ,APLL,APLC.APLR. 

2  APL . APCU . APLS . APCUS 

COMMON/D  I  MS/OPCUS  ,  HPCUS  .  DPCU  .  HPCU  .  DPLS  .  HPLS  .  DPL  .  HPL  , 

1  SD1.SD2.SD3 
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c 


c 

c 


c 

c 


c 


c 

c 


INITIALIZE  ALL  F'S  TO  0.0 
DO  5  L-1 ,N 
DO  10  M-1 ,N 
10  F(L.M)-0.0 
5  CONTINUE 

DETERMINE  VIEW  FACTOR  BETWEEN  PCUSR  AND  PLSL 

F( I  ST (10) . I  ST ( 12) )-DTOD(0.5>DPCUS,0 . 5. DPLS . SD 1 +S02+S03 ) 
F(IST(12), IST(i0))-F( I  ST (10) , I  ST ( 12) )*APCUSR/APLSL 


DETERMINE  VIEW  FACTORS  BETWEEN  PCU  AND  PCU  SHELL  SURFACES 
I  5-1  ST (5) 

I  7— I  ST ( 7 1 
19—1  ST  <  9 ) 

I  1 7— 1  ST (17) 

SD-(HPCUS-HPCU)»0.5 
FA-DTOD(0 . 5.DPCU , 0 . 5.DPCUS , SD) 

CALL  CTOR I NG(@ . 5.DPCU , HPCU , 0 . 5.DPCU . 0 . 5.DPCUS . SO . FB) 
FB^AMAXl (0.0.FB) 

F  '  ■  “  -  - 

F 
F 


117, I5J-APCUL/APCU*FA+APCUC/APCU*F8 
117, I9)-F( 117,15) 

117, I7)-2.0*APCUL/APCU*( 1 . 0-FA)+APCUC/APCU. ( 1 . 0-2.0.F8) 
15,1 1 7 )=APCU/ APCUSL*  F ( 117,15) 


’  I  9 ;  I 17)-F( 15, I  17) 

17,1 1 7 )—APCU/APCUSC*  F ( 117,17) 

ACCOUNT  FOR  PCU  BLOCKAGE 
B-(DPCUS— DPCU)»0.5 
D«B»(HPCUS/(HPCU+SD)-1 .0) 

D-AMIN1  (D,0.5»DPCU) 

DP=D+B 

CALL  FRINGS(0.5*DPCUS-DP,0.5*DPCUS.0.5*DPCU,0. 5.DPCUS , 


HPCUS , F ( 15,19)) 

19. 1  5 ) — F (  15,19) 

15. 1  7 )  — 1 . 0— F ( 15.1 1 7 )— F ( 15,19) 

19. 1  7 )-F(  15,  17) 

17, I5)-APCUSL/APCUSC«F( 15,17) 

17. 1  9 ) — F ( 17,15) 

17, I7)-1 ,0-2.0*F( 17, I5)-F( 17. I  17) 


DETERMINE  VIEW  FACTORS  BETWEEN  PL  AND  PL  SHELL  SURFACES 
I  5- 1  ST (5+6) 

I  7—1  ST (7+6) 

I  9—1  ST (9+6) 

1 1 7- 1  ST 7 17+1 ) 

SD-(HPLS-HPl)«0.5 

FA-DTOD ( 0.5.DPL, 0.5. DPLS, SD) 

CALL  CTOR I NG(0 . 5.DPL . HPL . 0 . 5.DPL , 0 . 5. DPLS , SD , FB) 
F8-AMAX1 (0.0, FB) 

F(  I  17  ,  I5)-APLL/APL*FA+APLC/APL.FB 
F( 117  .  I  9 )— F ( 117,15) 

F(  117  ,  I7)-2,0.APLL/APL»( 1 . 0-FA)+APLC/APL. ( 1 . 0-2.0. FB) 
F( 15,  I 17)-APL/APLSL.F( I  17,  15) 


F  (  I  9  .  I  1  7  )— F(  I  5  ,  117) 

F ( 1 7 , I 17)-APL/APLSC«F( 117,17) 

C  ACCOUNT  FOR  PL  BLOCKAGE 
B=(DPLS-OPL).0.5 
D=B* (HPLS/(HPL+SD)-1 .0) 

D-AMIN1 (D.0.5.DPL) 

DP-D+B 

CALL  FR  I  NGS(0 . 5* DPLS— DP  ,  0 . 5«DPLS  ,  0 . 5. DPL  ,  0 . 5* DPLS  , 

1  HPLS,F( 15, 19)) 

F ( 19. I  5 )— F ( 15.19) 

F( 15, I  7 >-1 ,0-F( 15, 1 17)-F( 15, 19) 

F( 19.  I  7 ) — F (  15.171 

F( 17, I5)-APLSL/APLSC«F( 15,17) 

F( 17, I  9)— F( 17,15) 

F( I  7, I  7 )  —  1  . 0-2,0. F( 17,1 5 )— F (17,117) 

C 

RETURN 

END 

C 

c  — *  *  *  * — 

c 

SUBROUT  t  NE  V  I EW2 ( N , NMR , NSR .  I  ST , RM , AM , RS . AS ) 

C  CALCULATE  VIEW  FACTORS  FOR  ALL  REMAINING  SURFACES 

DIMENSION  I  ST (20) , RM(NMR+1 ) , AM( NMR ) ,RS(NSR+1 ) , AS (NSR) 
COMMON/CV I EW/F (100,100) 

COMMON/ARE A/APCUSL , APCUSC , APCUSR , APCUL , APCUC , APCUR , 

1  APLSL ,APLSC,APLSR, APLL , APLC , APLR , 

2  APL , APCU , APLS , APCUS 

COMMON/D  I MS/DPCUS .HPCUS , DPCU . HPCU , DPLS , HPLS , DPL . HPL , 

1  SD1 .SD2.SD3 

C 

PI-3. 141592654 
C 

C  DETERMINE  VIEW  FACTORS  BETWEEN  MAIN  AND  SECONDARY  RADIATORS 
K— I  ST (2) 

DO  15  L-1 , NMR 
KK— I  ST (3) 

DO  20  M—1 .NSR 

CALL  FR I NGS (RM(L) ,RM(L+1 ) , RS (M) ,RS(M+1 ) , SD2 , FL ) 

FL-AMAX 1 ( 0 . 0 , FL ) 

F  (K ,  KK)«=FL 
C  REC I  PROS  I TY 

F(KK,K)-FL*AM(L)/AS(M) 

20  KK-KK+1 
15  K-K+1 
C 

C  DETERMINE  VIEW  FACTORS  BETWEEN  RIGHT  END  OF  PCU  SHELL  AND 
C  SECONDARY  RADIATOR 
KK= I  ST (3) 

DO  25  M-1 , NSR 

CALL  FR I NGS (0 .0.0. 5*DPCUS , RS (M) ,RS(M+1 ) , SD1+SD2 , FM) 
FM-AMAX1 (0.0.FM) 
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C  IF  SOI  IS  GREATER  THAN  ZERO,  MUST  ACCOUNT  FOR  BLOCKAGE  FROM 
C  MAIN  RADIATOR 
FB-0.0 

I F ( SD 1 . GT . 0 . 0 ) THEN 

B— SD1  /(SD1+SD2)*(RS(M+1  )-0.5*DPCUS)+0.5*DPCUS 
IF(B.GT.RM(1))THEN 

CALL  FR I NGS (6.0,0 . 5*DPCUS , RM( 1 ) ,0 , SD1 , F8) 

FB-AMAX1(0.0.FB) 

ELSE 
END  IF 
ELSE 
END  I  F 

F(  I  ST (10) ,KK)-AMAX1 (0.0.FM-FB) 

F  ( KK  ,  I  ST ( 10) )-FM»APCUSR/AS(M) 

25  KK-KK+1 
C 

C  DETERMINE  VIEW  FACTORS  BETWEEN  LEFT  END  OF  PL  SHELL  AND 
C  MAIN  RADIATOR 
KK-IST(2) 

DO  35  M-1 , NMR 

CALL  FR I NGS ( 0 , 0.0. 5*DPL$ , RM(M) ,RM(M+1 ) , SD2+SD3 , FM) 

FM-AMAX1 (0.0, FM) 

C  IF  SD3  IS  GREATER  THAN  ZERO,  MUST  ACCOUNT  FOR  BLOCKAGE  FROM 
C  SECONDARY  RADIATOR 
FB-0.0 

IF(SD3.GT.0.0)THEN 

B-=SD3/( SD2+SD3)  *  (RM(M+1  )-0 . 5*DPLS)+0 . 5«0PLS 
IF(B.GT.RS(1))THEN 

CALL  FRINGS(0.0,0.5.DPLS,RS(1) ,B,SD3,FB) 

FB-AMAX1(0.0,FB) 

ELSE 
END  IF 
ELSE 
END  IF 

F( I  ST ( 12) ,KK)-AMAX1 (0.0.FM-FB) 

F(KK,  I  ST ( 12) )-FM»APLSL/AM(M) 

35  KK-KK+1 
C 

C  DETERMINE  VIEW  FACTORS  BETWEEN  PCU  SHELL  CURVED  SURFACE 
C  AND  MAIN  RADIATOR 
K-IST(I) 

DO  45  L— 1 , NMR 

CALL  CTORING(0.5»DPCUS.HPCUS,RM(L).RM(L+1) ,SD1 , FL ) 

FL— AMAX 1 (0.0.FL) 

F ( I  ST ( 8 ) , K )— FL 
F ( K , IST(8))-FL.APCUSC/AM(L) 

45  K-K+1 
C 

C  DETERMINE  VIEW  FACTORS  BETWEEN  PCU  SHELL  AND  SECONDARY  RADIATOR 
C  IF  VIEW  IS  NOT  BLOCKED  BY  MAIN  RADIATOR 
RB=(HPCUS+SD 1 )/( HPCUS+SD 1 +SD2 ) 

K— I  ST ( 3 ) 


DO  50  L— 1 , NSR 
B-RB*RS(L+1 ) 

I F (B . LT . RM( 1 ) ) THEN 

CALL  CTORING(0.5*OPCUS,HPCUS,RS(L),RS(L+1), SD 1 +SD2 , FL ) 
FL«AMAX1(0.0,FL) 

F( IST(8) ,K)-FL 

F(K, IST(8))-FL*APCUSC/AS(L) 

ELSE 
END  IF 
50  K-K+1 
C 

C  DETERMINE  VIEW  FACTORS  BETWEEN  PL  SHELL  CURVED  SURFACE 
C  AND  SECONDARY  RADIATOR 
K—  I  ST  (4 ) 

DO  55  L— 1 , NSR 

CALL  CTORING(0. 5*DPLS , HPLS , RS ( L ) ,RS(L+1 ) .SD3.FL) 

FL-AMAX 1 (0.0. FL) 

F ( I  ST ( 1 4) , K )— FL 
F  (  K  ,  I  ST ( 1  4 ) )=FL*APL SC/AS ( L ) 

55  K-K+1 
C 

C  DETERMINE  VIEW  FACTORS  BETWEEN  PAYLOAD  SHELL  AND  MAIN  RADIATOR 
C  IF  VIEW  IS  NOT  BLOCKED  BY  SECONDARY  RADIATOR 
R8-(HPLS+SD3)/( HPLS+SD2+SD3 ) 

K— 1ST (2 ) 

DO  60  L— 1 , NMR 
B-RB*RM( L+1 ) 

IF(B.LT.RS(1))THEN 

CALL  CTOR  I  NG  ( 0 . 5*DPLS  .  HPLS  .  RM(  L)  ,  RM(  L+1  )  ,  SD2+SD3  ,  FL) 

FL— AMAX 1 ( 0 . 0 , FL ) 

F ( I  ST ( 1 4 ) , K )— FL 

F(K,  I  ST ( 1 4 ) )— FL*APLSC/AM( L ) 

ELSE 
END  IF 
60  K-K+1 
C 

C  IF  SD1  GREATER  THAN  ZERO,  DETERMINE  VIEW  FACTOR  BETWEEN 
C  PCU  SHELL  (RIGHT  END)  AND  MAIN  RADIATOR 
I F ( SD 1 . GT . 0 . 0 ) THEN 
K-IST(I) 

DO  65  L-1 , NMR 

CALL  FRINGS(0 . 0,0. 5*DPCUS,RM(L)  ,RM(L+1  )  ,  SD1  ,  FL  ) 

FL— AMAX  1  (  0 . 0  ,  FL  ) 

F ( I  ST ( 1 0) , K)— FL 

F(K, IST(10))-FL*APCUSR/AM(L) 

65  K-K+1 
ELSE 
END  I  F 
C 

C  IF  SD3  GREATER  THAN  ZERO,  DETERMINE  VIEW  FACTOR  BETWEEN 
C  PAYLOAD  SHELL  ( LETF  END)  AND  SECONDARY  RADIATOR 
IF(SD3.GT.0.0)THEN 
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no  non 


K-IST(4) 

CALL^FR  I  NGS(0 .0.0. 5>DPLS  ,RS  (L)  , RS ( L+1  )  .SD3.FL) 
FL-AMAX1 (0.0, FL) 

F ( 1 2) )— FL+APLSR/AS(L) 

75  K-K+1 
ELSE 
END  IF 

VIEW  FACTORS  BETWEEN  ALL  SURFACES  AND  SPACE 
MAIN  RADIATOR 
K-IST12) 

J-IST(l) 

DO  85  L-1  .  NMR 

SUM-0.0 

KK-IST(3) 

DO  90  M-1 , NSR 
SUM-SUM+F  ( K , KK ) 

90  SUM-SUM+F (K. IST( 12) )+F(K . IST( 14) ) 

^5:!!^11|)|:]1:f(j.ist(i0))-f(j.ist(8)) 
K-K+1 
85  J-J  +  1 

SECONDARY  RADIATOR 
K-IST13 
J-IST(4) 

DO  95  L-1 .NSR 
SUM-0 . 0 
KK— I  ST (2 ) 

DO  100  M— 1 , NMR 
SUM— SUM+F (K , KK ) 

100  SUM-SUM+F(K.IST(10))+F(K,IST(8)) 

IST(12))-F(J.  .ST(14)) 

K-K+1 
95  J-J+1 

PCU  SHELL  ,  „ 

F ( I  ST (6 ) ,  IST(19))-1 -0 
J— I  ST ( 3) 

SUM 1—0 . 0 

SUM2-0 . 0 

DO  105  L-1 .NSR 

SUM1— SUM1  +  F ( I  ST  f 1 0)  ,  J  ) 

SUM2— SUM2+F ( I  ST (8) . J) 

105  J-J+1 

K-IST(I) 

DO  1 10  M-1 .NMR 

SUM1— SUM1+F  C  ISTM  0)  ,K) 

SUM2-SUM2+F ( I  ST ( 8 ) ,K) 


F( l S  T ( 10) . IST(19))  — 1 . 0-SUM 1-F (IST(10),IST(12)) 
F( I  ST (8) .  I  ST (19) )-l  0-SUM2 
PL  SHELL 

F<  iItI^6) ’ IST( 0 

SUM1=02^ 

SUM2-00 

DO  115  L-1  .NMR 

SUM1— SUM1+F ( I ST( 12)  .  J) 

SUM2-SUM2+F ( I  ST ( 1  * ) , J ) 

115  J-J+1 

K—  I  ST(  4) 

DO  120  M-1 , NSR  ,  v  . 

SUM1— SUM1  +  F  ( IST( 12) ,K) 

SUM2-SUM2+F ( I  ST (14), K) 

120  ^ s4 (12).IST(19)}— 1  0-SUM1-F ( I  ST ( 1 2 ) . I  ST ( 10) ) 
F  ( 1  ST ( 1 4 ) , I  ST  (  1  9 )  )  — 1 . 0— SUM2 
SPACE 

F ( I  ST ( 1 9) .  IST(19))  =  1 .0 


SUBROUTINE  FASS I GN( I  FLAG .  I  ST , FP) 

ASSIGN  VIEW  FACTORS  FOR  PCU  AND  PL  TO  USE  IN 
SOLV I NG  COMB  I  NED  RAD  I  AT  I  ON/CONDUCT  I  ON 
COMMON/CV I EW/F (100.100) 

DIMENSION  FP(4 .4) . IST(20) 

I F ( I  FLAG . EO . 1 ) THEN 
FOR  PAYLOAD 
L-1  1 
M-1  8 
ELSE 
FOR  PCU 
L-5 
M-1  7 
END  IF 

K-IST(L) 

DO  5  1-1 . 4 
I F ( I  . EQ ■ 4 )K— I  ST (M) 

KK-IST(L) 

DO  10  J-1 ,4 

IF(J.E0.4)KK-IST(M) 

FP ( I .J)-F(K.KK) 

KK-KK+2 
10  CONTINUE 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


K-K+2 

5  CONTINUE 
RETURN 
END 


DETERM  ?NE°D  I  Sk'o^RADMJS^RI  TO  0  I! 5K  OF  RADIUS  R2 
WITH  H  SEPARATION  DISTANCE  (REF.  P . 826  #21  OF 
H-AMAX1 (H, 0.0001  ) 

CR1-R1/H 

CR2-R2/H 

DTOD-0.0 

I F(CR1 . LT . 1 ,0E-6)RETURN 
X-1 .0+(1 . 0+CR2*CR2)/CR1/CR1 
Y-CR2/CR1 

DTOD-0. 5* (X-(X«X-4 .0*Y«Y)*«0 . 5) 

RETURN 

END 


VIEW  FACTOR 
SIEGEL  AND  HOWELL 


) 


DETERM?NE°CYLIND^R1T026f^  VIEW  FACTOR  (USES  VIEW  FACTOR  ALGEBRA 
ALONG  WITH  #28 ,  P . 828  OF  SIEGEL  AND  HOWELL) 

R 1  1$  CYLINDER  RADIUS,  R2  IS  DISK  RADIUS 
CYLINDER  IS  ADJACENT  AND  PERPENDICULAR  TO  DISK 
Z  IS  LENGTH  OF  CYLINDER 
R-R2/R1 

IF(R.LT. 1 .0001 )THEN 

CTOD-0.0 

RETURN 

ELSE 

END  IF 

XL-Z/R1 

A-XL*XL+R*R-1 .0 

13  ( A+2 . 0)  **2-4 .0*R*R)**0. 5*  ACOS(B/R/A)+B*  AS  I  N(  1  .  0/R) 

1-1 .570796327*A 
C2-AC0S(B/A)-C1 *0 • 5/XL 
C3-1 .0/R-C2/R/3. 141592654 
F 1 2*=C3*R 

CTOD-0.5*(1 . 0-F12) 

RETURN 

END 


SUBROUTINE  FR I NGS(R1 I , RIO , R2 l ,R20,H,F) 
DETERMINE  VIEW  FACTORS  BETWEEN  RINGS  BASED 
VIEW  FACTORS  AND  VIEW  FACTOR  ALGEBRA 


ON  DISK  TO  DISK 


C 

C 

C 

c 

c 

C 

C 

c 

c 

c 

c 

c 


c 

c 

c 


c 


c 


c 


c 

c 

c 


F=(R10*R10* (DT0D(R10, 

1  R 1  I  * R 1  I  * (DT0D(R1 I , 

2  (R10*R10-R1 I *R1  I  ) 
RETURN 

END 


R20.H 

R20,H 


-DTOD ( R 1 0 , R2 I ,H 
-DTOD ( R 1  I  . R2  I  ,  H 


»)/ 


SUBROUT  I NE  CTOR I NG ( R 1  , Z 1 . R2 I  , R20.H , F) 
DETERMINE  VIEW  FACTORS  BETWEEN  CYLINDER 
R1  IS  CYLINDER  RADIUS.  Z1  IS  CYLINDER 
R2I  IS  RING  INNER  RADIUS.  R20  IS  RING 
H  IS  SEPARATION  DISTANCE 
Z1-AMAX1 ( Z 1 ,0.0001 ) 

EXTENDED  CYLINDER  TO  SMALL  RADIUS  DISK 
FF1-CT0D(R1 , R2 I ,Z1+H) 

CYLINDER  TO  SMALL  RADIUS  DISK 

FFF1-H/Z1*(FF1-CT0D(R1 .R2I . Z 1 ) )+FF1 
EXTENDED  CYLINDER  TO  LARGE  RADIUS  DISK 
FF2=CT0D(R1 .R20.Z1+H) 

CYLINDER  TO  LARGE  RADIUS  DISK 

FFF2-H/Z 1 . ( FF2-CT0D(R 1 , R20 , Z 1 ) )+FF2 
CYLINDER  TO  RING 
F-FFF2-FFF1 
RETURN 
END 


AND  RING 
LENGTH 

OUTER  RADIUS 


SUBROUTINE  OUT  1 (N) 

COMMON/EM  I  S/EM  I S( 100)  ,  , 

COMMON/HT/ONETPP (100) ,T(100) . R( 1 00) ,G( 1 00) 

WR [ TE ( 7  1 50 } 

150  FORMAt( :1 12X, ’NETWORK  SOLUTION  USING  GAUSS-S I EDEL ’ . 

1/  /.IX,  SURFACE  EMISSIVITY  TEMPERATURE  ONETPP  . 

2  6X ,  ' RAD  I  OS  I TY  I RR AO  I  AT  I  ON '/.22X,' (KELVIN)  ,4X,  (W/SO.M)  , 

3  6X, '(W/SO.M)' ,4X, ’ (W/SO.M)1 ,/) 


WRITE(7,160)(  I  ,EMIS(  I  )  .  T  (  I  )  .  ONETPP  Q)  ,  R  (  I  )  A  G(  O  .  I -1  .N) 
- 1 . -  —  "~^G.5,5X,F7.2.4X.1PE11 .4,2" 


160  FORMAT (4X, I2.2X.0PF9. 

RETURN 
END 


■2X.E1 i . 4 ,2X,E1 1 .4) 


SUBROUTINE  0UT2 ( 1ST) 

COMMON^AREA/APCUSU APCUSC , APCUSR , APCUL . APCUC . APCUR , 

1  / APLSL .APLSC.APLSR . APLL , APLC . APLR , 

2  APL.APCU.APLS.APCUS 
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i K ): ' ./. 

SHELL  - 


' ,F8.2,/. 


COMMON/OUT /PMR , PSR , PPCUS , PPCU , PPLS , PPL , PCUWH , PLWH , REFPOWE , 

1  PSR 1  , PSR2 , ARE AMR , AREASR , TAVGPCUS , T  AVGPLS , COP , 

2  T I SMR , TOSMR , T I SSR , TOSSR , HAPCU , HAPL , CPCU , CPL , 

3  THPOWER , THPNEW 
COMMON/HT /ONETPP ( 100) ,T (100) , R( 1 00) , G( 1 00 ) 

DIMENSION  I  ST (20) 

WR  I  TE  (6 , 1 00)  T  (  I  ST  (  1  7  )  )  ,  T  (  I  ST  (  6)  )  ,  T  (  I  ST  (8)  )  ,  T  (  I  ST  ( 1 0)  )  , TAVGPCUS 

100  FORMAT  (’1  1  , //, 4X ,  'PCU  TEMPERATURES  ' 

1 6X  ,  '  PCU  -  -  ,3X.F8.2,/,6X,  'LEFT  PCU 
26X , 1  CENTER  PCU  SHELL  -  ’,F8.2,/, 

36X,  'RIGHT  PCU  SHELL  -  \F8.2,/, 

46X, ’SHELL  AVERAGE  -  '  F8.2) 

WRITE (6, 105)T( I  ST (18)) ,T( I  ST (12)) ,T( I  ST (14)) ,T( I  ST (16)) .TAVGPLS 
105  FORMAT ( ’  ’ ,//,4X, ’PAYLOAD  TEMPERATURES  ( K )  :  ‘  ,  /  , 

16X, ’PAYLOAD  -  ’ , 3X , F8 . 2 ./ , 6X , ’ LEFT  PAYLOAD  SHELL  -  ',F8.2,/. 
26X, ’CENTER  PAYLOAD  SHELL  -  ’,F8.2,/. 

36X, ’RIGHT  PAYLOAD  SHELL  -  ',F8.2,/. 

46X, ’SHELL  AVERAGE  -  ’,F8.2) 

WRITE(6,107)hAPCU,HAPL 

107  FORMATC  ’ ,//.4X, ’TOTAL  HA  FOR  CONVECTION  (W/K):'./, 

16X,  ’PCU  ACTIVE  COOLING  -  ’  ,1  PE  12. 4./. 

26X, ’PAYLOAD  ACTIVE  COOLING  -  ’.E12.4) 

WCPCU— AMAX1  (CPCU  ,  XCPCU) 

WCPL-AMAX1 (CPL.XCPL) 

WRITE(6, 108)WCPCU,WCPL 

108  FORMATC  ’ ,//,4X, ’TOTAL  KA/X  FOR  CONDUCTION  (W/K):’,/. 

1 6X , ’ PCU  TO  SHELL  -  ’.1PE12.4,/. 

26X.  ’PAYLOAD  TO  SHELL  -  ’,E12.4) 

WRITE (6, 1 10)AREAMR. AREASR, APCUS.APCU.APLS.APL 
110  FORMATC  ’ ,//,4X, ’STRUCTURE  AREAS  (SO.M):  ' ,/, 

16X.-MAIN  RADIATOR  -  ’.1PE12.4,/, 

26X , ’SECONDARY  RADIATOR  -  ’.E12.4,/, 

36X , 'PCU  SHELL  -  ’ , E 1 2 . 4 ,/, 

46X, 'PCU  -  ’ ,E12.4,/, 

56X, 'PAYLOAD  SHELL  -  ’.E12.4./, 

66X , ’PAYLOAD  -  ’ .E12.4) 

WR I TE(6 , 1 1 5 ) PCUWH , PSR 1 , PPCUS 

115  FORMATC  ’ ,//.4X, ’WASTE  HEAT  GENERATED  IN  PCU  (W)  -  ’.1PE12.4. 
1/ , 6X , ’ PCU  WASTE  HEAT  REMOVED  ACTIVELY  (W)  -  ’.E12.4,/, 

26X ,  *  PCU  WASTE  HEAT  REMOVED  BY  CONDUCTION  TO  SHELL  (W)  -  ’  . 
3E12.4) 

WR 1 TE ( 6 . 120 ) PLWH .PSR2.PPLS 

120  FORMATC  ’ ,//,4X , ’WASTE  HEAT  GENERATED  IN  PAYLOAD  (W) 

1/.6X,  ’PAYLOAD  WASTE  HEAT  REMOVED  ACTIVELY  (W)  ■=  '.E12.4,/, 

26X,  ’PAYLOAD  WASTE  HEAT  REMOVED  BY  CONDUCTION  TO  SHELL  (W)  - 
3E12.4) 

WRITE (6, 125) PMR, PSR, PPCUS. PPCU. PPLS, PPL 
125  FORMATC  ’ ,//,4X, ’WASTE  HEAT  REMOVED  BY  RADIATION  (W) : ’  ./, 
16X,’MAIN  RADIATOR  -  ’.1PE12.4,/, 

26X , ’SECONDARY  RADIATOR  -  ’.E12.4,/, 

36X.’PCU  SHELL  -  ,E12.4,/,6X, ’ PCU  -  ’.E12.4,/, 

46X, ’PAYLOAD  SHELL  -  ’ , E 1 2 . 4 , / , 6X , ’ PAYLOAD  -  ’.E12.4) 


,  1  PE  1 2  . 


4  , 


POWMAX=AMAX 1 (THPOWER . THPNEW) 

WR I TE (6 . 1 30 ) REFPOWE , COP . POWMAX 

130  FORMATC  ’ ,//,4X , ’REFRIGERATOR  ELECTRIC  POWER  TO  COOL  PAYLOAD 
1  (W)  -  ’ ,1PE12. 4,/, 4X, ’REFRIGERATOR  COP  -  ’.E12.4./, 

24X, ’TOTAL  THERMAL  POWER  (W)  -  ’.E12.4,/. 

36X./.6X.25C  .  ’),//) 

C 

RETURN 

END 

C 

c  — ***♦ — 

c 

SUBROUTINE  I DOUT (N . NMR , NSR , IDSURF, IDCON) 

DIMENSION  IDSURF(n) , IDCON(N) 

WR I TE ( 7 , 1 00 ) 

100  FORMATC  ’ ,4X, ’SURFACE  IDENTIFICATION’,//, 

16X, 'IDSURF  -  0  -  SURFACE  WITH  KNOWN  TEMPERATURE',/, 

26X,’ IDSURF  -  1  -  SURFACE  WITH  KNOWN  HEAT  FLUX’./, 

36X,’ IDSURF  -  2  -  TWO-SIDED  SURFACE  WITH  T(SIDE  I)  - 

4  T ( S I DE  2)  AND ’ ,/,21X, 'Q(TOTAL)  •  QfSIDE  1)  +  0 ( S I DE  2)’./. 

56X, ’(IDCON  IDENTIFIES  CONNECTED  SURFACES,  0  FOR  NO  CONNECTION)’, 
0///.7X, ’SURFACE' ,8X, ' IDSURF' ,9X, ’ IDCON’ ,/) 

C 

DO  10  1-1 , N 

I  F ( I  . LE -  2*  NMR ) WR ITE(7, 105)1,  IDSURF(  I  )  ,  IDCON( I ) 

I  F(  I  .GT .2* NMR. AND.  I  . LE . 2*NMR+2»NSR) 

1 WRITE (7, 110) I . I DSURF ( I ) , IDCON( I ) 

I  F  (  I  . GT . 2»  NMR+2*  NSR )WR I TE ( 7 , 1 15) I  ,  I DSURF ( I ) . !DCON( I ) 

10  CONT INUE 

105  FORMAT?’  ’ ,7X, I3.12X, 13, 12X, I3.7X, ’MAIN  RADIATOR  NODE') 

110  FORMAT ( ’  ’ ,7X , 13, 12X,  13. 12X  .  13. 7X ,  ’SECONDARY  RADIATOR  NODE’) 

115  FORMAT ( ’  ’ ,7X, I3.12X, I3.12X, 13) 

RETURN 

END 

C 

c  — *  *  *  * — 

c 

SUBROUTINE  FOUT(N) 

COMMON/CV I EW/F (100,100) 

WR I TE (7 , 100) 

100  FORMAT( *  1 '  ,//,6X . ‘NON-ZERO  VIEW  FACTORS ‘ 

19X, ‘  I  *  ,6X , * J ‘ ,7X  ,  *F( I . J) ' ,6X, ‘F(J, 1 ) *  ,/) 

DO  10  1-1 , N 
DO  20  J-l ,N 

I F (ABS(F ( \ ,J)) . GT. 0.00005. OR. ABS(F(J,  l)) .GT. 0.00005) 

1WRITE(7, 105)1 ,J,F(I,J),F(J, I) 

20  CONTINUE 
10  CONTINUE 

105  FORMAT ( ’  ‘ ,6X, I3.4X, I3,6X,F7.4,5X,F7.4) 

RETURN 

END 
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—  *  *  *  * — 

SUBROUT  I NE  ROUT ( NMR . NSR , TNMR , TNSR . TSMR , TSSR , POAMR . POASR , 

1  RMR ,  RSR , AMR , ASR , RNMR , RNSR ) 

COMMON/OUT /PMR , PSR , PPCUS , PPCU , PPLS , PPL , PCUWH , PLWH , REFPOWE , 

1  PSR 1 , PSR2 , ARE AMR , AREASR , TAVGPCUS , TAVGPLS , COP , 

2  T  1  SMR  .  TOSMR  ,  T  I  SSR  ,  TOSSR  ,  HAPCU  ,  HAPL  ,  CPCU  ,  CPL  , 

3  THPOWER, THPNEW 

DIMENSION  TNMR(NMR) , TNSR ( NSR ) , T SMR (NMR ) .TSSR(NSR) 

DIMENSION  POAMR (NMR) .POASR (NSR) , RMR ( NMR+1 ) , RSR ( NSR+ 1 ) 
DIMENSION  AMR(NMR) .ASR(NSR) , RNMR (NMR ), RNSR ( NSR) 

WR I TE ( 6 , 1 00 ) RMR ( 1 ) , RMR ( NMR+ 1 ) . T I SMR , TOSMR 
100  FORMAT ( ’  ’ ,3X, ’MAIN-RADIATOR  NODE  OUTPUT',/, 

1 1 IX, ' INSIDE/OUTSIDE  RADII  (M)  -  ' , F8 . 3 , ' / ' , F8 . 3 , / , 

21 IX, ' INSIDE/OUTSIDE  TEMPERATURES  (K)  -  '  , F8 . 3 ,  ' / '  , F8 . 3 ,// . 
35X, 'NODE' ,3X, 'RADIUS  (M) ' . 3X , ' AREA  (SO .M) ' , 4X , ' TEMP  (K)', 
44X,'TEFF  (K) ' ,6X, 'POWER  (W)',/) 

DO  10  1-1, NMR 
P— POAMR ( I )*AMR( I ) 

10  WRITE (6, 115) I , RNMR ( I ) , AMR ( I ) , TNMR ( I ) ,TSMR( I ) , P 

WR I TE ( 6 , 1 05 ) RSR ( 1 ) , RSR( NSR+1 ) , T I SSR . TOSSR 
105  FORMAT ( ’  ',/,/,/, 3X, 'SECONDARY-RADIATOR  NODE  OUTPUT',/, 

1 1 1 X , ' INSIDE/OUTSIDE  RADII  (M)  -  ' , F8 . 3 , '/' , F8 . 3 ,/. 

211X, ' INSIDE/OUTSIDE  TEMPERATURES  (K)  -  ' , F8 . 3 , '/ ' , F8 . 3,//, 
35X, 'NODE' ,3X, 'RADIUS  (M) ' , 3X , ' AREA  ( SO . M) ' , 4X , ' TEMP  ( K ) ' , 
44X,'TEFF  (K) ' ,6X , 'POWER  (W)',/) 

DO  20  1-1 , NSR 
P-POASR( I ) *ASR( I ) 

20  WRITE (6. 1 15) I ,RNSR( I ) ,ASR( I ) ,TNSR( I ) , TSSR ( 1 ) ,P 

115  FORMAT ( '  ' ,4X, I3,4X,2(F8.3,5X),1X,2(F8.2,4X),1PE12.4) 

WR I TE (6 , 1 20 ) 

120  FORMAT ( '  -,//. 10X,20( •*-•') ,//) 

RETURN 
.  END 


C 

c  —  *  *  *  * — 
c 
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APPENDIX  C  -  SAMPLE  PROGRAM  OUTPUT 


MAIN-RADIATOR  NODE  OUTPUT 

INSIDE/OUTSIDE  RADII  (M)  =  6.000/  25.353 

INSIDE/OUTSIDE  TEMPERATURES  (K)  =  480.000/  985.000 


NODE  RADIUS  (M)  AREA  (SQ.M) 


TEMP  (K)  TEFF  (K) 


POWER  (W) 


1  17.233  1366.934 

2  23.244  381.674 

3  24.909  157.607 


549.20  372.66  9.9943E+06 
723.57  362.89  9.9986E+06 
893.99  360.90  9.9994E+06 


SECONDARY-RADIATOR  NODE  OUTPUT 

INSIDE/OUTSIDE  RADII  (M)  =  7.000/  44.572 

INSIDE/OUTSIDE  TEMPERATURES  (K)  =  470.000/  400.000 


NODE 

RADIUS  (M) 

AREA  (SQ.M) 

TEMP  (K) 

TEFF  (K) 

POWER  (W) 

1 

16 . 180 

1407.115 

460.63 

419.32 

2.0240E+06 

2 

26.728 

1433.203 

443 . 20 

396 . 58 

2.0268E+06 

3 

34.431 

1528.808 

425.73 

375.45 

2.0201E+06 

4 

41 . 242 

1718.179 

408.24 

356 . 91 

2.0196E+06 

*-**_**.**-**- 

A*-**-**-**-** 

-**.**-*•-** 

-**-**-**_** 

_**_**_* 
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pcu  TEMPERATURES  (K) : 

PCU  -  488 . 12 

LEFT  PCU  SHELL  -  459.17 

CENTER  PCU  SHELL  -  489.85 

RIGHT  PCU  SHELL  -  487.53 

SHELL  AVERAGE  -  485.75 


PAYLOAD  TEMPERATURES  (K)  : 

PAYLOAD  -  474.75 

LEFT  PAYLOAD  SHELL  -  508.12 

CENTER  PAYLOAD  SHELL  -  468.43 

RIGHT  PAYLOAD  SHELL  -  455.46 

SHELL  AVERAGE  -  471.80 


TOTAL  HA  FOR  CONVECTION  (W/K) : 

PCU  ACTIVE  COOLING  -  O.OOOOE+OO 
PAYLOAD  ACTIVE  COOLING  -  2.0378E+05 


TOTAL  KA/X  FOR  CONDUCTION  (W/K) : 
PCU  TO  SHELL  -  2.4514E+05 

PAYLOAD  TO  SHELL  -  4.4273E+05 


STRUCTURE  AREAS  (SQ.M) : 

MAIN  RADIATOR  -  1.9062E+03 

SECONDARY  RADIATOR  -  6. 087 3E +03 

PCU  SHELL  -  4.0212E+02 

PCU  -  2 . 4504E+02 

PAYLOAD  SHELL  -  6 . 2832E+02 

PAYLOAD  -  4.0212E+02 


WASTE  HEAT  GENERATED  IN  PCU  (W)  =  5.0000E+05 

PCU  WASTE  HEAT  REMOVED  ACTIVELY  (W)  =  >  O.OOOOE+OO 
PCU  WASTE  HEAT  REMOVED  BY  CONDUCTION  TO  SHELL  (W)  = 


4.9029E+05 


WASTE  HEAT  GENERATED  IN  PAYLOAD  (W)  =  9.0000E+06 

PAYLOAD  WASTE  HEAT  REMOVED  ACTIVELY  (W)  = 

PAYLOAD  WASTE  HEAT  REMOVED  BY  CONDUCTION  TO  SHELL  (W)  — 


8 . 8545E+05 


WASTE  HEAT  REMOVED  BY  RADIATION  (W)  : 
MAIN  RADIATOR  -  3.0000E+07 

SECONDARY  RADIATOR  -  8.1000E+06 

PCU  SHELL  -  4 . 9029E+05 

PCU  -  1 . 4547E+04 

PAYLOAD  SHELL  -  8.8545E+05 

PAYLOAD  -  1 . 4547E+04 
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SURFACE  IDENTIFICATION 

IDSURF  =  O  -  SURFACE  WITH  KNOWN  TEMPERATURE 
IDSURF  =  1  -  SURFACE  WITH  KNOWN  HEAT  FLUX 

IDSURF  =  2  -  TWO-SIDED  SURFACE  WITH  T(SIDE  1)  =  TfSIDE  2)  AND 

Q (TOTAL)  =  Q(SIDE  1)  +  Q(SIDE  2) 

(IDCON  IDENTIFIES  CONNECTED  SURFACES;  O  FOR  NO  CONNECTION, 
NEGATIVE  FOR  SECOND  SIDE  OF  TWO-SIDED  SURFACE) 

SURFACE  IDSURF  IDCON 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 
29 


O 

O 

O 

O 

O 

0 

0 

O 

O 

0 

O 

0 

O 

O 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

2 

1 

1 

O 


4 

5 

6 
-1 
-2 
-3 
11 
12 

13 

14 
-7 
-8 
-9 

-10 

16 

-15 

18 

-17 

20 

-19 

22 

-21 

24 

-23 

26 

-25 

0 

O 

O 


MAIN  RADIATOR  NODE 
MAIN  RADIATOR  NODE 
MAIN  RADIATOR  NODE 
MAIN  RADIATOR  NODE 
MAIN  RADIATOR  NODE 
MAIN  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 
SECONDARY  RADIATOR  NODE 


PCU  SHELL  LEFT  END  (INSIDE) 
PCU  SHELL  LEFT  END  (OUTSIDE) 
PCU  SHELL  CENTER  (INSIDE) 

PCU  SHELL  CENTER  (OUTSIDE) 
PCU  SHELL  RIGHT  END  (INSIDE) 
PCU  SHELL  RIGHT  END  (OUTSIDE) 
PL  SHELL  LEFT  END  (INSIDE) 

PL  SHELL  LEFT  END  (OUTSIDE) 
PL  SHELL  CENTER  (INSIDE) 


PL  SHELL  CENTER 


PL  SHELL  RIGHT  END 
PL  SHELL  RIGHT  END 
PCU 


(INSIDE) 

(OUTSIDE) 


PAYLOAD  (PL) 
SPACE 
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NON- ZERO  VIEW  FACTORS 


I 

J 

F(I,J) 

F(J,I) 

1 

18 

0.0543 

0.2462 

1 

20 

0.0016 

0.0438 

1 

29 

0.9441 

0.0000 

2 

18 

0.0156 

0.0197 

2 

20 

0.0001 

0.0005 

2 

29 

0.9843 

0.0000 

3 

18 

0.0127 

0.0066 

3 

20 

0.0000 

0.0001 

3 

29 

0.9873 

0.0000 

4 

7 

0.2478 

0.2408 

4 

8 

0.1650 

0.1574 

4 

9 

0.1155 

0.1033 

4 

10 

0.0862 

0.0686 

4 

22 

0.0167 

0.2899 

4 

29 

0.3687 

0.0000 

5 

7 

0.1921 

0.0521 

5 

8 

0.1613 

0.0429 

5 

9 

0.1268 

0.0317 

5 

10 

0.0998 

0.0222 

5 

22 

0.0108 

0.0525 

5 

29 

0.4093 

0.0000 

6 

7 

0.1774 

0.0199 

6 

8 

0.1576 

0.0173 

6 

9 

0.1287 

0.0133 

6 

io 

0.1035 

0.0095 

6 

22 

0.0097 

0.0194 

6 

29 

0.4232 

0.0000 

7 

20 

0.0104 

0.2906 

7 

29 

0.6769 

0.0000 

8 

20 

0.0055 

0.1575 

8 

29 

0.7768 

0.0000 

9 

20 

0.0034 

0.1023 

9 

29 

0.8484 

0.0000 

10 

20 

0.0022 

0.0747 

10 

29 

0.8976 

0.0000 

11 

22 

0.0020 

0.0364 

11 

24 

0.0743 

0.2218 

11 

29 

0.9237 

0.0000 

12 

22 

0.0001 

0.0010 

12 

24 

0.0193 

0.0587 

12 

29 

0.9807 

0.0000 

13 

22 

0.0000 

0.0004 

13 

24 

0.0096 

0.0311 

13 

29 

0.9904 

0.0000 

14 

22 

0.0000 

0.0002 

14 

24 

0.0057 

0.0208 

14 

29 

0.9943 

0.0000 

15 

17 

0.3321 

0.0553 

15 

19 

0.0366 

0.0366 

15 

27 

0.6313 

0.1295 

16 

29 

1.0000 

0.0000 

17 

17 

0.2872 

0.2872 

17 

19 

0.0553 

0.3321 

17 

27 

0.6021 

0.7410 
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18 

29 

0.7274 

0.0000 

19 

27 

0.6313 

0.1295 

20 

22 

0.0235 

0.0150 

20 

29 

0.3071 

0.0000 

21 

23 

0.3226 

0.0538 

21 

25 

0.0298 

0.0298 

21 

28 

0.6476 

0.1265 

22 

29 

0.5852 

0.0000 

23 

23 

0.2550 

0.2550 

23 

25 

0.0538 

0.3226 

23 

28 

0.6375 

0.7470 

24 

29 

0.6677 

0.0000 

25 

28 

0.6476 

0.1265 

26 

29 

1.0000 

0.0000 

29 

29 

1.0000 

1.0000 

NETWORK 

SOLUTION  USING  GAUSS-SIEDEL 

EMISSIVITY 

TEMPERATURE 

QNETPP 

RADI OS I TY 

(KELVIN) 

(W/SQ.M) 

(W/SQ.M) 

0.90000 

549 . 20 

4.0939E+03 

4.7027E+03 

0.90000 

723.57 

1 . 3534E+04 

1 . 4035E+04 

0.90000 

893.99 

3 . 2146E+04 

3 . 2639E+04 

0.90000 

549 . 20 

3 . 2170E+03 

4.8001E+03 

0.90000 

723.57 

1 . 2662E+04 

1.4132E+04 

0.90000 

893.99 

3 . 1299E+04 

3 . 2733E+04 

0.90000 

460.63 

-2 . 9950E+02 

2.5854E+03 

0.90000 

443 . 20 

-1.0456E+02 

2 . 1988E+03 

0.90000 

425.73 

7.6644E+01 

1 . 8538E+03 

0.90000 

408.24 

1 . 8230E+02 

1 . 5544E+03 

0.90000 

460.63 

1.7373E+03 

2 . 3591E+03 

0.90000 

443 . 20 

1 . 5186E+03 

2.0185E+03 

0.90000 

425.73 

1 . 2447E+03 

1 . 7241E+03 

0.90000 

408.24 

9 . 9310E+02 

1.4643E+03 

0.90000 

459 . 18 

-6 . 3371E+02 

2 . 5906E+03 

0.90000 

459 . 17 

1 . 8547E+03 

2 . 3139E+03 

0.90000 

489.88 

6 . 9416E+01 

3 . 2572E+03 

0.90000 

489 . 85 

1 . 1506E+03 

3 . 1362E+03 

0.90000 

487.53 

6.6968E-01 

3. 2026E+03 

0.90000 

487.53 

1 . 2187E+03 

3.0673E+03 

0.90000 

508.13 

8 . 6508E+02 

3 . 6832E+03 

0.90000 

508.12 

5 . 4606E+02 

3 . 7184E+03 

0.90000 

468 . 47 

-1 . 1749E+02 

2 . 7435E+03 

0.90000 

468.43 

1.5276E+03 

2 . 5598E+03 

0.90000 

455 . 46 

-3.7299E+02 

2.4809E+03 

0.90000 

455 . 46 

1.7822E+03 

2 . 2414E+03 

0.80000 

488.12 

3 . 9644E+01 

3 . 2082E+03 

0.80000 

474.75 

3 . 6175E+01 

2.8707E+03 

1.00000 

300.00 

0 .  OOOOE+OO 

4.5919E+02 

IRRADIATION 

(W/SQ.M) 

6.0883E+02 
5.0108E+02 
4 . 9330E+02 
1.5831E+03 
1.4695E+03 
1.4347E+03 
2.8849E+03 
2 . 3034E+03 
1.7772E+03 
1 . 3721E+03 
6 . 2184E+02 
4 . 9990E+02 
4.7936E+02 
4.7121E+02 
3 . 2243E+03 
4.5919E+02 
3 . 1878E+03 
1 . 9855E+03 
3 . 2019E+03 
1 . 8486E+03 
2.8181E+03 
3 . 1724E+03 
2 . 8610E+03 
1.0323E+03 
2.8539E+03 
4.5919E+02 
3 . 1638E+03 
2.8292E+03 
4.5919E+02 
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3315  Old  Forest  Road 
P.O.  Box  10935 
Lynchburg,  VA  24506-0935 
Attn:  B.  J.  Short 

Battelle  Pacific  Northwest  Lab. 
P.  O.  Box  999 
Richland,  WA  99352 
Attn:  J.  O.  Barner 

Battelle  Pacific  Northwest  Lab. 
P.  O.  Box  999 
Richland,  WA  99352 
Attn:  L.  Schmid 

Battelle  Pacific  Northwest  Lab. 
P.  O.  BOX  999 
Richland,  WA  99352 
Attn:  E.  P.  Coomes 


Battelle  Pacific  Northwest  Lab 
Battelle  Boulevard 
Richland,  WA  99352 
Attn:  B.  M.  Johnson 

Battelle  Pacific  Northwest  Lab 
P.  O.  Box  999 
Richland,  WA  99352 
Attn:  W.  J.  Krotiuk 

Battelle  Pacific  Northwest  Lab 
P.  O.  Box  999 
Richland,  WA  99352 
Attn:  R.  D.  Widrig 

\ 

Boeing  Company 
P.O.  Box  3999 
MS  8K-30 

Seattle,  WA  98124-2499 
Attn:  A.  Sutey 

Boeing  Company 
Boeing  Aerospace  System 
P.O.  Box  3707 
Seattle,  WA  98124 
Attn:  K.  Kennerud 

Brookhaven  National  Laboratory 
P.O.  Box  155 
Upton,  NY  11973 
Attn:  T.  Bowden 

Brookhaven  National  Laboratory 
P.O.  Box  155 
Upton,  NY  11973 
Attn:  H.  Ludewig 

Brookhaven  National  Laboratory 
P.O.  Box  155 
Upton,  NY  11973 
Attn:  W.  Y.  Kato 

Brookhaven  National  Laboratory 

P.O.  Box  155 

Bldg.  701,  Level  143 

Upton,  NY  11973 

Attn:  J.  Powell 

California  Inst,  of  Technology 
Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 
Attn:  V.  C.  Truscello 


Dist-3 


California  Inst,  of  Technology 
Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 
Attn:  P.  Bankston 

California  Inst,  of  Technology 
Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 
Attn:  E.  P.  Framan 

California  Inst,  of  Technology 
Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 
Attn:  L.  Isenberg 

California  Inst,  of  Technology 
Jet  Propulsion  Laboratory 
4800  Oak  Grove  Drive 
Pasadena,  CA  91109 
Attn:  J.  Mondt 

nippi 

1400  Wilson  Blvd. 

Arlington,  VA  22209 
Attn:  P.  Kemmey 

DCS CON  Consulting 
4265  Drake  Court 
Livermore,  CA  94550 
Attn:  D.  C.  Sewell 

Defense  Nuclear  Agency 
6801  Telegraph  Road 
Alexandria,  VA  22310-3398 
Attn:  J.  Farber/RAEV 

DNA/RAEV 

6801  Telegraph  Road 
Alexandria,  VA  22310-3398 
Attn:  J.  Foster 

EG&G  Idaho,  Inc./INEL 
p.O.  Box  1625 
Idaho  Falls,  ID  83415 
Attn:  R.  Rice 


EG&G  Idaho,  Inc./INEL 
P.O.  BOX  1625 
Idaho  Falls,  ID  83415 
Attn:  J.  Dearien 

EG&G  Idaho,  Inc./INEL 
P.O.  Box  1625 
Idaho  Falls,  ID  83415 
Attn:  M.  L.  Stanley 

EG&G  Idaho,  Inc./INEL 
P.O.  Box  1625 
Idaho  Falls,  ID  83415 
Attn:  R.  D.  Struthers 

EG&G  Idaho,  Inc./INEL 
P.O.  Box  1625 
Idaho  Falls,  ID  83415 
Attn:  J.  F.  Whitbeck 

EG&G  Idaho,  Inc./INEL 
P.O.  Box  1625 
Idaho  Falls,  ID  83415 
Attn:  P.  W.  Dickson 

EG&G  Idaho,  Inc./INEL 
P.O.  Box  1625 
Idaho  Falls,  ID  83415 
Attn:  J.  W.  Henscheid 

Ford  Aerospace  Corporation 
Aeronutronic  Div. 

Ford  Road,  P.O.  Box  A 
Newport  Beach,  CA  92658-9983 
Attn:  V.  Pizzuro 

GA  Technologies 
P.O.  Box  85608 
San  Diego,  CA  92138 
Attn:  H.  J.  Snyder 

GA  Technologies 
P.O.  Box  85608 
San  Diego,  CA  92138 
Attn:  C.  Fisher 

GA  Technologies 
P.O.  Box  85608 
San  Diego,  CA  92138 
Attn:  R.  Dahlberg 
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Garrett  Fluid  Systems  Co. 
P.O.  Box  5217 
Phoenix,  AZ 
Attn:  Robert  Boyle 

General  Electric  Company 
P.  O.  Box  8555 
Astro  Systems 
Philadelphia,  PA  19101 
Attn:  R.  J.  Katucki 

General  Electric  Corp/NSTO 
310  DeGuigne  Drive 
Sunnyvale,  CA  9048 
Attn:  C.  Cowan 


Grumman  Aerospace  Corporation 
M/S  B20-05 
Bethpage,  NY  11714 
Attn:  J.  Belisle 

Hanford  Engineering  Dev.  Lab 
Post  Office  Box  1970 
Richland,  WA  99352 
Attn:  D.  S.  Dutt 

House  of  Representatives  Staff 
Space  and  Technology  Committee 
2320  Rayburn  Building 
Washington,  DC  20515 
Attn:  Tom  Weimer 


General  Electric  NSTO 
310  DeGuigne  Drive 
Sunnyvale,  CA  90486 
Attn:  E.  E.  Gerrels 

General  Electric  NSTO 
310  DeGuigne  Drive 
Sunnyvale,  CA  90486 
Attn:  H.  S.  Bailey 

General  Electric-SCO 
P.  O.  Box  8555 
Astro  Systems 
Philadelphia,  PA  19101 
Attn:  J.  Chan 

General  Electric 
P.  O.  Box  8555 
Bldg.  100,  Rm  M2412 
Astro  Systems 
Philadelphia,  PA  19101 
Attn:  J.  Hnat 


Idaho  National  Engineering 
Laboratory 
P.  O.  Box  1625 
Idaho  Fall,  ID  83414 
Attn:  W.  H.  Roack 

Innovative  Nuclear  Space  Pwr.Inst 
202  NSC 

University  of  Florida 
Gainesville,  FL  32611 
Attn:  N.  J.  Diaz 

International  Energy  Assoc.  Ltd. 
1717  Louisiana  NE 
Suite  202 

Albuquerque,  NM  87110 
Attn:  G.  B.  Varnado 

Lawrence  Livermore  National  Lab. 
P.  O.  Box  808 
Livermore,  CA  94550 
Attn:  Lynn  Cleland,  MS  L-144 


General  Electric 
P.  O.  Box  8555 
Astro  Systems 
Philadelphia,  PA  19101 
Attn:  R.  D.  Casagrande 

General  Electric 
P.  O.  Box  8555 
Astro  Systems 
Philadelphia,  PA  19101 
Attn:  W.  Chiu 


Lawrence  Livermore  National  Lab. 
P.  O.  Box  808 
Livermore,  CA  94550 
Attn:  C.  E.  Walter,  MS  L-144 

Los  Alamos  National  Laboratory 
P.  O.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  T.  Trapp,  MS-E561 
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Los  Alamos  National  Laboratory 
P.  O.  Box  1663 
Los  Alamos ,  NM  87545 
Attn:  R.  Hardie,  MS-F611 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  C.  Bell,  MS  A145 

Los  Alamos  National  Laboratory 
P.  O.  BOX  1663 
Los  Alamos,  NM  87545 
Attn:  R.  J.  LeClaire 

Los  Alamos  National  Laboratory 
P.  O.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  S.  Jackson,  MS-F611 

Los  Alamos  National  Laboratory 
P.  O.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  J.  Metzger 

Los  Alamos  National  Laboratory 

P.  0.  Box  1663 

Los  Alamos,  NM  87545 

Attn:  C.  W.  Watson,  MS-F607 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  Don  Reid,  MS-H811 

Los  Alamos  National  Laboratory 
P.  O.  BOX  1663 
Los  Alamos,  NM  87545 
Attn:  R.  Bohl,  MS-K551 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  D.  R.  Bennett 

Los  Alamos  National  Laboratory 
P.  O.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  W.  L.  Kirk 


Los  Alamos  National  Laboratory 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  M.  Merrigan 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  T.  P.  Suchocki 

Los  Alamos  National  Laboratory 
P.  0.  Box  1663 
Los  Alamos,  NM  87545 
Attn:  L.  H.  Sullivan 

Martin  Marietta  Corp. 

P.  0.  Box  179 
Denver,  CO  80201 
Attn:  R.  Giellis 

Martin  Marietta  Corp. 

P.  0.  Box  179 
Denver,  CO  80201 
Attn:  R.  Zercher 
MSL8060 

Massachusetts  Institute  of 
Technology 
1328  Albany  Street 
Cabridge,  MA  02139 
Attn:  J.  A.  Bernard 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  Barbara  McKissock,  MS  301-5 


NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  H.  Bloomfield,  MS  301-5 


NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  A.  Juhasz,  MS  301-5 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  J.  Smith,  MS  301-5 
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NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  C.  Purvis,  MS  203-1 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  D.  Bents,  MS  301-5 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  I.  Myers,  MS  301-2 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  G.  Schwarze,  MS  301-2 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  J.  Sovie,  MS  301-5 

NASA  Lewis  Research  Center 
21000  Brookpark  Road 
Cleveland,  OH  44135 
Attn:  Barbara  Jones,  MS  301-5 

National  Research  Council 
Energy  Engineering  Board 
Commission  on  Engineering 
and  Technical  Systems 
2101  Constitution  Avenue 
Washington,  DC  20418 
Attn:  R.  Cohen 

Naval  Research  Laboratory 
Washington,  DC  20375-5000 
Attn:  R.  L.  Eilbert 

Naval  Research  Laboratory 
Washington,  DC  20375-5000 
Attn:  I.  M.  Vitkovitsky 

Naval  Space  Command 
Dahlgren,  VA  22448 
Attn:  Commander  R.  Nosco 


Naval  Space  Command 
N5 

Dahlgren,  VA  22448 
Attn:  Maj.  J.  Wiley 

Naval  Space  Command 
Dahlgren,  VA  22448 
Attn:  Mr.  B.  Meyers 

Naval  Surface  Weapons  Center 
Dahlgren,  VA  22448-5000 
Attn:  R.  Gripshoven-F12 

Naval  Surface  Weapons  Center 
Dahlgren,  VA  22448-5000 
Attn:  R.  Dewitt-F12 

Naval  Surface  Weapons  Center 
White  Oak  Laboratory 
Silver  Springs,  MD  20903-500 
MC  R-42 

Attn:  B.  Maccabee 

Nichols  Research  Corp. 

2340  Alamo  Street,  SE 
Suite  105 

Albuquerque,  NM  87106 
Attn:  R.  Weed 

Oak  Ridge  National  Laboratory 

P.  O.  Box  Y 

Bldg.  9201-3,  MS-7 

Oak  Ridge,  TN  37831 

Attn:  J.  p.  Nichols 

Oak  Ridge  National  Laboratory 

P*  O.  Box  Y 

Bldg.  9201-3,  MS-7 

Oak  Ridge,  TN  37831 

Attn:  D.  Bartine 

Oak  Ridge  National  Laboratory 
P.  O.  Box  X 
Oak  Ridge,  TN  37831 
Attn:  H.  W.  Hoffman 

Oak  Ridge  National  Laboratory 

P.  O.  Box  y 

Bldg.  9201-3,  MS-7 

Oak  Ridge,  TN  37831 

Attn:  R.  H.  Cooper,  Jr. 
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Oak  Ridge  National  Laboratory 

p.  O.  Box  Y 

Bldg.  9201-3,  MS-7 

Oak  Ridge,  TN  37831 

Attn:  J.  C.  Moyers 

Oak  Ridge  National  Laboratory 
P.  O.  Box  Y 
Oak  Ridge,  TN  37831 
Attn:  M.  Olszewski 


SDI  Organization 
The  Pentagon 

Washington,  DC  20301-7100 
Attn:  R.  Wiley 

SDI/SLKT 
The  Pentagon 
1717  H.  St.  NW 
Washington,  D.  C.  20301 
Attn:  C.  Northrup 


Oak  Ridge  National  Laboratory 

P.  O.  Box  Y 

Bldg.  9201-3,  MS-7 

Oak  Ridge,  TN  37831 

Attn:  M.  Siman-Tov 

Oak  Ridge  National  Laboratory 

p.  O.  Box  Y 

Bldg.  9201-3,  MS-7 

Oak  Ridge,  TN  37831 

Attn:  F.  W.  Wiffen 

RADC/OCTP 
Griffiss  AFB 
New  York  13441 
Attn:  R.  Gray 

Riverside  Research  Institute 
1701  No.  Ft.  Meyers  Drive 
Suite  700 

Arlington,  VA  22209 
Attn:  J.  Feig 

Science  Applications,  Inc. 

505  Marquette  Avenue  NW 
Albuquerque,  NM  87102 
Attn:  D.  Buden 

Science  &  Engineering  Associates 
6301  Indian  School  Road,  NE 
Albuquerque,  NM  87110 
Attn:  G.  L.  Zigler 


SDIO/DE 

Washington,  DC  20301-7100 
Attn:  Dr.  J.  Hammond 

SDIO/IST 

Washington,  DC  20301-7100 
Attn:  Dr.  L.  Cavery 

SDIO/KE 
The  Pentagon 

Washington,  DC  20301-7100 
Attn:  Col.  R.  Ross 

SDIO/KE 
The  Pentagon 

Washington,  DC  20301-7100 
Attn:  Maj.  R.  X.  Lenard 

SDIO/SATKA 

Washington,  DC  20301-7100 
Attn:  Col.  Garry  Schnelzer 

SDIO/SY 

Washington,  DC  20301-7100 
Attn:  Dr.  C.  Sharn 

SDIO/SY 

Washington,  DC  20301-7100 
Attn:  Col.  J.  Schofield 

SDIO/SY 

Washington,  DC  20301-7100 
Attn:  Col.  J.  Graham 


SDI  Organization 


The  Pentagon 
Washington,  DC 
Attn:  R.  Verga 


20301-7100 


SDIO/SY 

Washington,  DC  20301-7100 
Attn:  Capt.  J.  Doegan 
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Space  Power,  Inc. 

253  Humbolt  Court 
Sunnyvale,  CA  94089 
Attn:  J.  R.  Wetch 

State  University  of  New  York 
at  Buffalo 

Dept,  of  Elec.  Engineering 
312  Bonner  Avenue 
Buffalo,  NY  14260 
Attn:  Jim  Sargeant 

TRW 

One  Space  Park 
Redondo  Beach,  CA  90278 
Attn:  R.  Hammel 

TRW 

One  Space  Park 
Redondo  Beach,  CA  90278 
Attn:  T.  Fitzgerald 

TRW 

One  Space  Park 
Redondo  Beach,  CA  90278 
Attn:  C.  Garner 

TRW-ATD 
One  Space  Park 
Redondo  Beach,  CA  90278 
Attn:  B.  Glasgow 

TRW 

One  Space  Park 
Redondo  Beach,  CA  90278 
Attn:  A.  D.  Schoenfeld 

Teledyne  Brown  Engineering 
Cummings  Research  Park 
Huntsville,  AL  35807 
Attn:  Dan  DeLong 

Texas  A&M  University 
Nuclear  Engineering  Dept. 
College  Station,  TX  77843-3133 
Attn:  F.  Best 

Texas  Tech,  university 
Dept,  of  Electrical  Engr. 
Lubbock,  TX  79409 
Attn:  Dr.  W.  Portnoy 


U.  S.  Army  ARDC 
Building  329 
Picatinny  Arsenal 
New  Jersey  87806-5000 
Attn:  SMCAR-SSA-E 

U.  S.  Army  Belvoir  RDE  Center 
Fort  Belvoir,  VA  22060-5606 
Attn:  Dr.  L.  Amstutz-STRABE-FGE 

U.  S.  Army  Lab.  Com. 

SLKET/ML 

Pulse  Power  Technology  Branch 
Fort  Monmouth,  NJ  07703-5000 
Attn:  S.  Levy 

U.  S.  Army  Lab.  Com. 

SLKET/ML 

Pulse  Power  Technology  Branch 
Fort  Monmouth,  NJ  07703-5000 
Attn:  N.  Wilson 

U.  S.  Army  Strategic  Defense  Com. 
106  Wynn  Drive 
Huntsville,  AL  35807 
Attn:  C.  Cooper 

U.  S.  Army  Strategic  Defense  Com. 
106  Wynn  Drive 
Huntsville,  AL  35807 
Attn:  G.  Edlin 

U.  S.  Army  Strategic  Defense  Com. 
106  Wynn  Drive 
Huntsville,  AL  35807 
Attn:  R.  Hall 

U.  S.  Army  Strategic  Defense  Com. 
106  Wynn  Drive 
Huntsville,  AL  35807 
Attn:  E.  L.  Wilkinson 

U.  S.  Army  Strategic  Defense  Com. 
106  Wynn  Drive 
Huntsville,  AL  35807 
Attn:  D.  Bouska 

U.  S.  Army  Strategic  Defense  Com. 
106  Wynn  Drive 
Huntsville,  AL  35807 
Attn:  W.  Sullivan 
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D.  S.  Army  Strategic  Defense  Com. 
106  Wynn  Drive 
Huntsville,  AL  35807 
Attn:  F.  King 

D.  S.  Department  of  Energy 
Chicago  Operations  Office 
9800  S.  Cass  Avenue 
Argonne,  IL  60439 
Attn:  J.  L.  Hooper 


U .  S .  Department  of  Energy 
Pittsburgh  Energy  Tech.  Center 
P.0.  Box  18288 
Pittsburgh,  PA  15236 
Attn:  G.  Staats  (PM-20) 

U.  S.  Department  of  Energy 
NE-54 

Washington,  DC  20545 
Attn:  I.  Helms 


D.  S.  Department  of  Energy 

NE-52 

GTN 

Germantown,  HD  20545 
Attn:  J.  Warren 

U.  S.  Department  of  Energy 

NE-54 

F415/GTN 

Germantown,  HD  20545 
Attn:  E.  Wahlquist 

U.  S.  Department  of  Energy 
NE-52 1 

Germantown,  MD  20874 
Attn:  D.  Bennett 

U.  S.  Department  of  Energy 
San  Francisco  Operations  Office 
1333  Broadway  Ave, 

Oakland,  CA  94612 
Attn:  J.  K.  Hartman 

U.  S.  Department  of  Energy 
SAN  -  ACR  Division 
1333  Broadway 
Oakland,  CA  94612 
Attn:  J.  Krupa 

D.  S.  Department  of  Energy 
SAN  -  ACR  Division 
1333  Broadway 
Oakland,  CA  94612 
Attn:  W.  Lambert 

U.  S.  Department  of  Energy 
SAN  -  ACR  Division 
1333  Broadway 
Oakland,  CA  94612 
Attn:  J.  Zielinski 


U.  S.  Department  of  Energy 
Oak  Ridge  Operations  Office 
P.0.  Box  E 
Oak  Ridge,  TN  37830 
Attn:  E.  E.  Hoffman 

U.  S.  Department  of  Energy 
Washington,  DC  20545 
Attn:  s.  j.  Lanes 

U.  S.  Department  of  Enerav 
MA  206 

Washington,  DC  20545 
Attn:  J.  p.  Lee 

U.  S.  Department  of  Energy 
San  Francisco  Operations  Office 
1333  Broadway  Avenue 
Oakland,  CA  94612 
Attn:  S.  L.  Samuelson 

U.  S.  Department  of  Enerav 

ALO/ETD 

P.O.  Box  5400 

Albuquerque,  New  Mexico  87115 
Attn:  R.  Holton 

U.  S.  Department  of  Enerav 

ALO/ETD 

P.O.  Box  5400 

Albuquerque,  New  Mexico  87115 
Attn:  c.  Quinn 

U.  S.  Department  of  Energy/Idaho 
785  DOE  Place 
Idaho  Falls,  ID  83402 
Attn:  p.  j.  Dirkmaat 
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United  Technologies 
International  Fuel  Cells 
195  Governor's  Highway 
South  Windsor,  CT  06074 
Attn:  D.  McVay 

United  Technologies 
International  Fuel  Cells 
195  Governor's  Highway 
South  Windsor,  CT  06074 
Attn:  J.  L.  Preston,  Jr. 

United  Technologies 
International  Fuel  Cells 
195  Governor's  Highway 
South  Windsor,  CT  06074 
Attn:  J.  C.  Trocciola 

University  of  Missouri  -  Rolla 
220  Engineering  Research  Lab 
Rolla,  MO  65401-0249 
Attn:  A.  S.  Kumar 

University  of  New  Mexico 
Chemical  and  Nuclear  Eng, 
Department 

Albuquerque,  NM  87131 
Attn:  M.  El-Genk 

University  of  Wisconsin 
Fussion  Technology  Institute 
1500  Johnson  Drive 
Madison,  WI  53706-1687 
Attn:  Gerald  Kukinski 

W.  J.  Schafer  Associates 
1901  No.  Ft.  Myers  Drive 
Suite  800 

Arlington,  VA  22209 
Attn:  P.  Mace 

W.  J.  Schafer  Associates 
1901  No.  Ft.  Myers  Drive 
Suite  800 

Arlington,  VA  22209 
Attn:  S.  Bassett 

W.  J.  Schafer  Associates 
1901  No.  Ft.  Myers  Drive 
Suite  800 

Arlington,  VA  22209 
Attn:  M.  Nikolich 


W.  J.  Schafer  Associates 
1901  No.  Ft.  Myers  Drive 
Suite  800 

Arlington,  VA  22209 
Attn:  J.  Crissey 

W.  J.  Schafer  Associates 
2000  Randolph  Road,  SE 
#205 

Albuquerque,  NM  87106 
Attn:  D.  C.  Straw 

W.  J.  Schafer  Associates 
1901  No.  Ft.  Myers  Drive 
Suite  800 

Arlington,  VA  22209 
Attn:  A.  K.  Hyder 

Westinghouse  Electric 
P.  O.  Box  158 
Madison,  PA  15663-0158 
Attn:  J.  Chi 

Westinghouse 

Advanced  Energy  Systems  Division 
Manager,  Space  &  Defense  Program 
Route  70,  Madison  Exit 
Madison,  PA  15663 
Attn:  J.  F.  Wett 

Westinghouse 

Advanced  Energy  Systems  Division 
P.O.  Box  158 
Madison,  PA  15663 
Attn:  Dr.  J.  W.  H.  Chi 

Westinghouse  R&D 
1310  Beulah  Road 
Bldg.  501-3Y56 
Pittsburgh,  PA  15235 
Attn:  J.  R.  Repp 

Westinghouse  R&D 
1310  Beulah  Road 
Bldg.  501-3Y56 
Pittsburgh,  PA  15235 
Attn:  L.  Long 


Dist-11 


Westinghouse  R&D 

6400 

D.  McCloskey 

1210 

Beulah  Road 

6410 

N.  Ortiz 

Bldg. 

501-3Y56 

6420 

J.  Walker 

Pittsburgh.  PA  15235 

6421 

P.  Pickard 

re 

Attn: 

Owen  Taylor 

6422 

J.  Brockman 

6425 

W.  Camp 

Westinghouse  Advanced  Energy 

6431 

J.  Philbin 

Systems  Division 

6440 

D.  Dahlgren 

P.  O. 

Box  158 

6450 

T.  Schmidt 

Madison.  PA  15663 

6500 

A.  W.  Snyder 

Attn: 

G.  Farbman 

6510 

W.  Gauster 

6511 

L.  Cropp 

Westinghouse  Hanford  Co. 

6511 

M.  Edenburn 

P.  0. 

Box  1970 

6511 

D.  Gallup 

Richland.  WA  99352 

6511 

S.  Hudson 

Attn: 

D.  S.  Dutt 

6511 

A.  Marshall 

6511 

W.  McCulloch 

Westinghouse  Hanford  Co. 

6511 

R.  Pepping 

P.  0. 

BOX  1970 

6511 

F.  Thome 

Richland.  WA  99352 

6512 

D.  Ericson  (5) 

Attn: 

B.  J.  Makenas 

6512 

M.  Chu 

6512 

V.  Dandini 

6512 

D.  Dobranich  (10) 

1140 

P.  Peercy 

6512 

P.  McDaniel 

1200 

J.  P.  Van  Devender 

6512 

F.  Wyant 

1240 

K.  Prestwich 

8024 

P.  W.  Dean 

1248 

M.  Buttram 

8400 

R.  Wayne 

1270 

R.  Miller 

9000 

R.  Hagengruber 

1271 

M.  Clauser 

9010 

W.  C.  Hines 

1512 

D.  Rader 

9012 

J.  Keizur 

1800 

R.  Schwoebel 

9012 

L.  Connell 

1810 

G.  Kepler 

9012 

R.  Zazworski 

1830 

M.  Davis 

9100 

R.  Clem 

1832 

W.  Jones 

9110 

P.  Stokes 

1832 

R.  Salzbrenner 

9140 

D.  Rigali 

1840 

R.  Eagan 

2110 

R.  Bair 

2120 

W.  Dawes.  Jr. 

2140 

C.  Gibbon 

2150 

E.  Graham.  Jr. 

2560 

J.  Cutchen 

3141 

S.  A.  Landenberger  (5) 

3151 

W.  L.  Garner  (3) 

3154- 

1  C.H.  Dalin.  for  DOE/OSTI 

(28) 

* 
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